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Abstract.  We  study  the  approximation  properties  of  a  harmonic  function 
u  £  k  >  0,  on  a  relatively  compact  subset  A  of  Q,  using  the  Gen¬ 

eralized  Finite  Element  Method  (GFEM).  If  £1  =  0,  for  a  smooth,  bounded 
domain  O,  we  obtain  that  the  GFEM-approximation  ug  £  S  of  u  satisfies 
||rr. — A  C/G  ||«|  -h  .  where  h  is  the  typical  size  of  the  “elements” 
defining  the  GFEM— space  S  and  7  >  0  is  such  that  the  local  approximation 
spaces  contain  all  polynomials  of  degree  k  +  7.  The  main  technical  ingredi¬ 
ent  is  an  extension  of  the  classical  super-approximation  results  of  Nitsche  and 
Schatz  [20,  21].  In  addition  to  the  usual  “energy”  Sobolev  spaces  H 1(0),  we 
need  also  the  duals  of  the  Sobolev  spaces  Hm(0),  m  £  Z_|_. 
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Introduction 

Let  us  consider  the  Neumann  problem 

(  A u  =  0  on  O, 

}  dvu  =  g  G  H'-'V'HdO)  on  dO, 


where  O  is  a  smooth,  bounded  open  subset  of  R",  dO  is  the  boundary  of  O,  and  du 
is  the  directional  derivative  in  the  direction  of  the  outer  unit  normal  v  to  dO.  In 
this  paper,  we  are  interested  mainly  in  the  case  r  <  1,  r  £  Z,  and  we  are  looking 
for  the  approximation  properties  of  the  solution  u  £  Hr(0)  on  suitable  subsets  of 
O. 
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For  v  G  Hr{0),  r  G  R,  r  >  3/2,  the  boundary  values  (or  traces)  v\ go  and 
dvv\ go  are  defined  classically,  because  the  restriction  to  the  boundary  extends  by 
continuity  to  maps  Hr(0)  3!)-*  v\ oo  €  Hr~1/2(dO)  and  Hr(0 )  3»-t  d„v\ oo  € 
Hr-V2(dO),  see  [11,  31]  for  example.  For  r  <  3/2,  this  is  no  longer  true  in  general, 
but  for  v  —  u,  we  can  take  advantage  of  the  fact  that  u  satisfies  an  elliptic  equation, 
so  it  is  still  possible  to  define  dvu  G  Hr~3/2(dO)  [15]  (see  also  [18,  28]).  We  can 
assume,  without  loss  of  generality,  that  O  is  connected,  that  is,  that  O  is  a  domain. 

It  is  then  known  [15,  26,  28]  that  a  solution  u  of  Equation  (1)  exists  for  any  g 
such  that  ( g ,  l)eo  =  0  and  that  this  solution  satisfies 

(2)  IM|jd-(0)  <  0\\g\\Hr- 3/2(9e)), 

with  a  constant  C  that  may  depend  on  r,  but  is  independent  of  g.  (Here  (g,  l)eo  is 
the  value  of  the  distribution  g  on  the  function  constant  equal  to  1,  thus  (g,  l)eo  = 
Sdo  9(x)dS(x)  if  g  is  a  function.)  In  particular,  to  solve  the  boundary  value  problem 
(1),  it  is  enough  to  do  that  for  gn  G  C°°(dO),  gn  — >  g  g  Hr~3^2(dO).  The  estimate 
of  Equation  (2)  will  be  discussed  in  detail  in  [6],  where  more  references  will  be  given. 
In  that  paper,  the  case  A u  ^  0  will  also  be  considered.  In  this  paper,  however,  we 
avoid  all  together  the  issue  of  defining  the  traces  u\qo  by  introducing  in  Section  4 
a  variational  formulation  of  the  boundary  value  problem  (1).  This  also  leads  to  a 
quick  derivation  of  Equation  (2).  See  also  [18]. 

The  consideration  of  the  case  when  g  g  Hr~3/2{0)1  r  <  3/2,  when  g  is  a 
distribution  rather  than  a  function,  is  important  in  order  to  be  able  to  handle  the 
case  of  “concentrated  moments”  and  “concentrated  loads,”  see  for  example  [27]. 
The  concentrated  loads  and  moments  are  distributions  concentrated  at  one  point, 
so  they  are  obtained  by  taking  derivatives  of  ^-distributions.  See  Example  1.4  and 
below  for  the  definition  of  the  delta  distributions  and  of  their  derivatives.  Most 
of  our  results  work  without  any  reference  to  boundary  conditions.  However,  the 
definition  of  the  numerical  (or  discrete  solution)  is  more  difficult  for  the  Dirichlet 
problem,  so  we  do  not  address  this  problem  explicitly  in  this  paper. 

Let  O  C  Rn  denote  a  bounded,  connected,  open  subset  of  R"  (i.e.,  will  be  a 
bounded  domain) .  We  do  not  assume  that  H  is  smooth,  unless  explicitly  mentioned 
(if  O  is  assumed  to  be  smooth,  then  we  shall  use  the  notation  O  instead  of  Q). 
Recall  that  A  <g  B  means  that  A  is  bounded,  is  contained  in  the  interior  of  B  and 
dA  and  dB  are  disjoint  (i.e.,  A  is  a  relatively  compact  subset  of  B).  If  A  <s  fl  is  an 
open  subset,  then  the  solution  u  of  Equation  (1)  will  be  smooth  on  A  for  any  r  and 

(3)  IM|iT™(A)  <  C,||u||Jyr(fl), 

with  a  constant  C  that  depends  on  A,  H,  r,  and  m,  but  is  independent  of  u  satis¬ 
fying  A u  =  0.  An  important  problem,  with  potential  practical  applications,  is  to 
approximate  on  A  the  solution  u  of  Equation  (1). 

In  this  paper,  we  prove  several  results  on  the  approximation  of  the  solution  u  on 
subsets  A  <g  O  using  the  Generalized  Finite  Element  Method.  Let  Sv  C  Hm(0), 
m  >  k  +  1,  be  a  sequence  of  Generalized  Finite  Element  Spaces  associated  to  a 
sequence  =  {wj,  <j>j,  \E/-, v  G  Z+,  of  GFEM-data  with  typical  size  hv 
satisfying  the  assumptions  of  Subsection  2.3  (so,  in  particular,  hv  — »  0).  Then 
Theorem  5.8  gives  that  the  sequence  of  GFEM-approximations  uu  :=  usv  £  Sv  of 
the  solution  u  of  the  boundary  value  problem  (1)  satisfies 

(4)  ll'ii  -  Uv II 1(A)  <  CKl\\u\\Hi-k{0), 
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provided  that  our  local  approximation  spaces  contain  all  polynomials  of  degree 

k  +  7,  7  >  0. 

The  idea  of  considering  partitions  of  unity  and  of  the  Generalized  Finite  Element 
Method  was  introduced  by  Babuska,  Caloz,  and  Osborn  [3].  It  was  further  devel¬ 
oped  in  [2,  5,  16],  and  [4].  The  Generalized  Finite  Element  Method  is  used  today 
in  Engineering  under  various  names,  such  as:  the  method  of  “clouds,”  the  method 
of  “finite  spheres,”  the  “X-finite  element  method,”  and  others.  The  Generalized 
Finite  Element  Method  is  a  generalization  of  the  mesh  free  (meshless)  methods 
which  use  as  a  paradigm  the  idea  of  partition  of  unity  introduced  in  [3,  5,  14]  and 
[16].  See  [2]  and  [14]  for  further  references. 

We  stress  that  our  results  require  not  just  the  energy  Sobolev  space  H1,  but  also 
negative  order  Sobolev  spaces  H~l ,  defined  in  this  paper  as  the  duals  of  Hl ,  l  £  Z+. 
One  of  the  main  reasons  for  the  need  to  consider  the  negative  order  Sobolev  spaces 
is  that  the  solution  u  is  in  H1~k(fl),  and  not  in  H1~k( ft),  in  general.  Moreover, 
even  if  we  approximate  the  boundary  data  g  and  the  solution  u  with  functions  in 
H1 ,  then  it  will  still  be  important  to  use  the  norm  on  a  negative  order  Sobolev 
space  in  the  estimate  of  the  error. 

Here  is  now  a  brief  description  of  the  contents  of  the  paper.  We  continue  to 
assume  that  S!  is  a  bounded  domain,  but  we  do  not  assume  that  fi  is  smooth,  except 
when  explicitly  mentioned.  A  domain  that  is  assumed  to  be  smooth  will  be  usually 
denoted  by  O .  In  Section  1,  we  set  up  the  notation  and  we  establish  our  conventions 
on  Sobolev  spaces.  Section  2  contains  a  quick  review  of  the  necessary  definitions 
involving  the  Generalized  Finite  Element  Method  (GFEM).  Our  main  results  are 
statements  on  a  sequence  Sv  of  GFEM-spaces  satisfying  the  assumptions  of  Sub¬ 
section  2.3.  Our  assumptions  are  formulated  in  terms  of  four  general  conditions 
(Conditions  A (h),  B,  C,  and  D),  formulated  in  Subsection  2.2.  The  spaces  Sv  will 
contain  the  sequence  uu  £  Sv  of  approximations  of  the  solution  u  to  our  bound¬ 
ary  value  problem  (Equation  (1)).  The  following  section,  Section  3,  contains  the 
calculations  necessary  to  establish  our  interior  estimates  for  the  sequence  u „  £  S„, 
v  £  Z+.  Our  approach  follows,  to  a  certain  extend,  that  in  the  article  of  Nitsche 
and  Schatz  [21],  relying  also  on  Wahlbin’s  survey  article  [33].  The  main  differences 
between  the  approach  in  Section  3  of  our  paper  and  the  approach  in  [21,  33]  to  inte¬ 
rior  estimates  are  due  mostly  to  the  fact  that  several  assumptions  from  those  papers 
are  not  fully  satisfied  in  our  approach.  As  in  those  articles,  the  main  step  is  a  super¬ 
approximation  property,  Proposition  3.7.  The  proof  in  [21,  33]  cannot  be  used  to 
obtain  Proposition  3.7  because  the  property  “ daw  =  0  if  |a|  is  large,”  is  not  satisfied 
in  general  for  w  £  Sv.  For  the  results  of  Sections  4  and  5,  we  assume  that  Q  =  O  is 
smooth.  In  Section  4,  we  introduce  a  weak  formulation  of  the  Neumann  problem  (1) 
and  the  Galerkin  approximation  (or  GFEM-approximation)  u$  £  S  of  the  solution 
u.  This  is  based  on  an  extension  B  of  the  form  B(w,v)  :=  f0Vu(x )  •  \/v(x)dx 
to  the  case  when  v  £  Hk~l(JD)  is  arbitrary  and  u  £  Hl~k(<D)  can  be  written  as 
u  =  Ui  +  U2,  where  Au±  =  0  in  distribution  sense  and  w2  G  i71(0).  We  also  es¬ 
tablish  a  well  posedness  result  for  (1)  by  establishing  that  the  extension  B  satisfies 
the  Babuska-Brezzi  condition.  In  Section  5  we  exploit  the  definition  and  properties 
of  the  Galerkin  approximations  uv  :=  us„  £  Sv  and  of  the  form  B.  Several  esti¬ 
mates  for  u  £  H1~k(0)  and  its  approximations  u„  are  established  in  this  section, 
including  the  main  theorem,  Theorem  5.8  (whose  main  conclusion  was  summarized 
in  Equation  (4)  above).  The  last  section,  Section  6  contains,  in  particular,  a  proof 
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that,  for  a  domain  Cl  with  piecewise  C[ -boundary,  we  can  construct  a  family  of 
partitions  of  unity  Sv  with  typical  size  of  supports  hv  — >  0  that  satisfies  the  as¬ 
sumptions  of  our  main  results  ( i.e .,  the  assumptions  of  Subsection  2.3)  for  a  fixed 
choice  of  the  structural  constants  (i.e.,  of  A,  Cj,  a ,  k,  A,  and  to).  For  this  con¬ 
struction,  we  assume  that  the  local  approximation  spaces  are  'Fj  =  Q\ ,  the  space 
of  polynomials  of  degree  at  most  A.  By  contrast,  it  is  not  possible  to  find  a  family 
of  partitions  of  unity  as  above  for  domains  with  cusps,  see  Remark  6.7.  For  suitable 
g  G  H~1^2~k(dO),  we  plan  to  perform  some  concrete  numerical  simulations  in  a 
future  paper  [7],  where  we  shall  also  give  more  examples  of  domains  and  partitions 
of  unity  satisfying  he  Assumption  (iv)  of  Subsection  2.3. 

We  shall  write  x  :=  y  if  x  is  defined  by  y.  By  C  we  shall  denote  a  constant 
that  may  depend  only  on  the  dimension  n  (O  C  R").  By  A ,  Cj,  a,  k,  A,  and 
m,  we  shall  denote  the  so  called  “structural  constants”  introduced  in  Conditions 
A-D  (Subsection  2.2).  The  structural  constants  will  remain  fixed  throughout  our 
discussion.  By  contrast,  C  will  denote  a  generic  constant  that  may  depend  only 
on  the  structural  constants  (and,  occasionally,  on  subsets  A,  A',...  C  Cl,  when 
explicitly  mentioned) . 

The  second  named  author  thanks  G.  Grubb,  A.  Schatz,  and  L.  Wahlbin  for 
some  useful  references.  A.  Schatz  has  also  made  some  useful  comments  on  an 
earlier  version  of  the  manuscript,  for  which  we  are  grateful. 

1.  Preliminaries 

We  begin  by  fixing  the  notation  and  terminology.  We  denote  by  R  the  set  of 
real  numbers  and  by  C  :=  {a  +  bi,a,b  £  M}  the  set  of  complex  numbers.  Also, 
N  =  {1,  2, . . .}  and  Z+  =  {0}  UN.  By  L2(Cl)  we  shall  denote  the  space  of  bounded, 
measurable  functions  /  :  — >  C  with  norm  II/II^q)  :=  fn  \f(x)\2dx.  (As  usual, 
we  identify  two  functions  that  coincide  outside  a  set  of  Lebesgue  measure  zero.) 

We  shall  use  the  multi-index  notation  for  partial  derivatives.  Namely,  a  = 
(<*1, . . . ,  an)  €  Z"  will  denote  a  multi-index  and  dau  :=  2  . . .  <9“" u,  as  usual. 

Also,  |a|  :=  «i  +  02  +  •  •  •  +  ,  is  the  order  of  the  partial  derivative  da. 

From  now  on  we  shall  assume  that  Cl  is  a  domain,  that  is,  an  open,  connected 
subset  o/Rn. 

1.1.  Sobolev  spaces.  We  now  introduce  the  integral  Sobolev  spaces  on  Cl.  This 
is  enough  for  our  purposes.  Let  s  £  Z+.  Then  Hs(Cl)  is  the  space  of  functions 
/  €  L2(Cl)  such  that 

(5)  11/112=  E  II^/IIIro)  <  °o- 

|a|<s 

The  space  Hq(CI)  is  defined  as  the  closure  of  C£°(Cl)  in  Hs (Cl). 

We  define  the  negative  order  Sobolev  spaces  by  duality.  Namely,  H~s(Cl)  := 
Hs(Cl)*,  the  dual  of  Hs(Cl),  s  £  N.  We  shall  denote  by  (w,<j>)  =  w(fi)  the  value 
of  the  linear  functional  w  £  H~s(Cl )  on  </>  £  Hs(Cl).  Then  the  norm  on  H~s(Cl)  is 
given  by 

IM|g-(n)=sup  ,  0  ^(j>£Hs(Ci). 

Our  definition  of  negative  order  Sobolev  spaces  by  duality  follows  [10,  22,  25], 
for  example.  Note,  however,  that  the  negative  order  Sobolev  spaces  are  often 
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also  defined  by  restriction  from  R",  as  in  [11,  15,  31],  for  example.  The  space 
of  restrictions  to  of  distributions  in  H~s(M.n)  is  the  dual  of  Hq(Q),  and  will  be 
denoted  HqS( fi).  The  spaces  HqS( Q)  =  s  >  0,  will  also  be  used  below. 

When  SI  =  R™,  these  two  approaches  yield  the  same  spaces,  but  for  general  0  they 
may  lead  to  different  “negative  order”  Sobolev  spaces.  The  spaces  Hs(M.n),s  £  R, 
can  also  be  defined  using  the  Fourier  transform. 

1.2.  Distributions.  Since  Cf°(Q)  C  Hs( SI)  for  any  s  €  Z+,  we  obtain  that  every 

w  £  s  £  Z+,  defines  a  distribution  on  SI.  For  s  <  0,  the  spaces  Hs( SI) 

consist,  in  general,  of  distributions  and  not  of  functions.  We  shall  not  use  any 
nontrivial  results  on  distributions,  but  we  shall  use  the  terminology  related  to  dis¬ 
tributions,  so  we  now  review  a  few  needed  definitions.  Let  B r(0)  denote  the  open 
ball  of  radius  R  centered  at  the  origin.  Also,  let  C£°(R")  be  the  set  of  infinitely 
differentiable,  complex  valued  functions  that  vanish  outside  a  ball  Br( 0),  for  some 
large  R  >  0.  The  elements  of  this  space  are  sometimes  called  test  functions.  A 
linear  map  u  :  C“(R”)  — >  C  is  called  a  distribution  on  R”  [12,  13,  31]  if,  for  any 
R  >  0,  there  exists  m  £  Z+  and  C  >  0  such  that 

(6)  \u(<j>)\  <  C  ^  ||9“(/>||l2,  if  (f  £  C“(R11)  and  <j>  =  0  outside  Br( 0). 

|  ck|  <ra 

This  definition  does  not  exclude  the  case  when  larger  and  larger  values  of  m  and  C 
have  to  be  chosen  as  R  — »  oo,  and  in  fact  this  situation  actually  occurs  in  specific 
examples.  The  set  of  distributions  on  R"  will  be  denoted  2?'(R"). 

We  now  fix  more  notation  and  terminology.  If  /  is  a  function,  then  the  closure  of 
the  set  {/  yf  0}  is  called  the  support  of  f  and  will  be  denoted  supp(/).  Therefore, 
any  <j)  £  C^° (fl)  has  compact  support.  We  shall  also  write  (u,  <j>)  :=  u{(j>)  for  the  value 
of  the  distribution  u  on  the  function  </>  £  C£°(Rn).  The  support  of  a  distribution  u 
is  the  smallest  closed  set  F  such  that  ( u ,  (f>)  =  0  for  any  </>  £  C“(Rn  \  F). 

Here  are  some  examples  of  distributions  and  constructions  based  on  distributions 
that  are  useful  below. 

Example  1.1.  If  /  is  a  measurable  function  on  R™  that  is  integrable  on  any  closed 
ball  in  R”  ( i.e it  is  locally  integrable,  or  /  £  L11oc(R")),  then  we  can  define 

(7)  (/,</>)  :=  f  f{x)cf(x)dx, 

J  R" 

for  any  </>  £  C£°(R”).  Thus  any  /  £  L11oc(Rn)  defines  a  distribution  on  R”,  that  is, 

AoC(R”)  C  2?'(R”). 

The  derivatives  of  distributions  are  defined  by  duality. 

Example  1.2.  The  derivatives  dau  of  a  distribution  u  are  defined  by 

<■ dau,cf )  :=  (-1)IQI (u,dact>). 

1.3.  Sobolev  spaces  on  the  boundary.  We  shall  also  need  the  definition  of 
the  spaces  Hrn+1'2{dO)  for  O  a  smooth,  bounded  domain  and  m  £  Z+.  Then 
Hrn+1/2(ydO)  consists  of  the  restrictions  to  dO  of  the  functions  in  Hm+1{0)  with 
norm 


(8) 


IMI  Hm+1/2(aO)  —  llwllfl-m+rfo) 
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where  w  G  Hm+1(0 )  is  the  unique  solution  of  Aw  =  0  and  w\go  =  v  (m  G  Z+).  An 
equivalent  norm  is  given  by  inf  ||u||_f/m+i(Ci),  where  u  G  Hm+l(0)  satisfies  u  =  v  on 
DO.  If  r  =  m+ 1/2,  m  G  Z+,  the  space  H~r(dO)  is  defined  as  the  dual  of  Hr(dfl). 
We  shall  denote  by  (u,  4>)dO  =  v{4>)i  the  value  of  v  G  H~r(dO)  on  </>  G  Hr{dO). 

Example  1.3.  If  n  =  2,  then  each  connected  component  of  90  is  diffeomorphic  to 
S'1,  the  unit  circle.  It  is  therefore  enough  to  define  I?r(S1)  in  this  case,  which  has  a 
more  concrete  description.  Let  /  G  C°°(S1).  Then,  up  to  a  multiplicative  constant 

OO 

(9)  \\f\\H^)=^  E  l/(n)|2(l+M)2r, 

n=— oo 

where  2tt f(n)  =  f(e2m0)e~2mnffdd,  n  G  Z,  are  the  Fourier  coefficients  of  /.  Then 
Hr(S1)  is  the  closure  of  C°°(S1)  in  the  norm  ||/||_Hr(si). 

The  following  example  is  relevant  for  the  discussion  of  concentrated  loads  and 
moments. 

Example  1.4.  The  Dirac  measure  (or  distribution)  concentrated  at  a  G  80  is  the 
distribution  Sa  defined  by 

(SaA)dO  :=^(o),  4>  G  C°°(dO). 

An  explicit  calculation  shows  that  5a  G  H~(n~~1^2-e(dO)  and  the  resulting  norms 
behave  as  ||5a||ij-("-i)/2-«  ~  e-1/2  — *  oo  as  e  — »  0. 

In  Section  4,  it  will  be  convenient  to  use  another  definition  of  fractional  Sobolev 
spaces  on  dO ,  which  yields  a  more  suitable  form  of  the  inner  product.  Namely,  let 
A  be  the  Laplace  operator  on  dO.  Recall  that  A  is  defined  by 

—  (Av,v)  =  f  |Vr)|2d5(a;). 

JdO 

Then  let  Uj  G  L2(dO),  j  G  Z+,  be  an  orthonormal  basis  of  L2(dO)  consisting 
of  eigenfunctions  of  —A,  that  is,  —A  Uj  =  XjUj.  We  can  assume  that  A  j  <  \j+i 
for  all  j  G  Z_|_.  In  particular,  it  follows  that  Ao  =  0  and  that  uq  is  constant  on 
each  connected  component  of  dO.  We  then  define  Hs'2{dO)  to  be  the  set  of  those 
functions  u  =  aj  uj  e  L2{dO)  for  which 

IMI^(OO)  :=  £(!  +  **)  \aj I  <  °°- 

j 

If  also  v  =  bjuj  G  Hs(dO),  then  the  inner  product  in  Hs(dO)  is  given  by 

(u,v)ffs(dO)  =  E(X  +  -\?)SaA- 
j 

The  advantage  of  this  formula  is  that  it  is  immediately  seen  that 

(10)  (u,  l)ff.(ao)  =  /  u(x)dS(x), 

JdO 

for  all  s  >  0.  For  s  <  0,  the  same  conclusion  follows  by  duality. 
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2.  The  Generalized  Finite  Element  Method 

We  now  recall  a  few  basic  facts  about  the  Generalized  Finite  Element  Method 
[2,  5,  16].  This  method  is  quite  convenient  when  one  needs  finite  element  spaces 
with  high  regularity.  Most  of  the  results  of  this  section  work  for  a  general  bounded 
open  set  f 2,  except  the  application  to  the  boundary  value  problems,  Section  4,  in 
which  case  we  shall  need  to  assume  that  f 2  has  a  smooth  boundary. 


2.1.  Basic  facts.  Let  k  €  Z+.  We  shall  denote  as  usual 

Mw*.°°(n)  :=  max  ||<9Qu||Loo(n),  ||u||Wk,oo(n)  :=  max  ||dau||Loo(n), 

|  a  | = Ac  |  o:  |  <  A; 

Wfc’°°(M)  :=  {u,||u||w*.»(n)  <  oo},  and  || Vw|| Wk,oo{n)  :=  T,j\\dM\w^(n)-  In 
particular,  \u\wo,<»(n)  =  ||u||w°>~(Q)  =  |M|l°°(Q)- 

We  shall  need  the  following  slight  generalization  of  a  definition  from  [5,  16]: 

Definition  2.1.  Let  O  C  R"  be  an  open  set  and  be  an  open  cover  of  O 

such  that  any  ie!l  belongs  to  at  most  k  of  the  sets  ojj.  Also,  let  {<pj}  be  a  partition 
of  unity  consisting  of  W/m’°°(fl)  functions  and  subordinated  to  the  covering  {u^-} 
(i.e.,  supp^j  C  cJj).  If 

(11)  l|9Q^||Loo(0)  <  C'fc/(diamwi)fc,  k  =  \a\  <  m, 

for  any  j  =  1, . . . ,  N,  then  {4>j}  is  called  a  ( k ,  Co,  C i, . . . ,  Cm)  partition  of  unity. 

Assume  also  that  we  are  given  linear  subspaces  Tj  G  Hm(uij),  j  =  1,2, ...  ,7V. 
The  spaces  will  be  called  local  approximation  spaces  and  will  be  used  to  define 
the  space 

N 

(12)  S  =  Sqfem  ■=  |  vj  €  C  Hm(Lkjjj), 

i=i 

which  will  be  called  the  GFEM-space.  The  set  {uJj,  <pj,  Tj}  will  be  called  the  set  of 
data  defining  the  GFEM-space  S. 

A  basic  approximation  property  of  the  GFEM-spaces  is  the  following  Theorem 
from  [5]. 


Theorem  2.2  (Babuska-Melenk).  We  shall  use  the  notations  and  definitions  of 
Definition  2.1  and  after.  Let  {f>j}  be  a  (k,Cq,Ci)  partition  of  unity.  Also,  let 
Vj  €  Tj  C  uap  :=  'Yhj4>jvj  £  >5,  and  dj  =  diam  u)j,  the  diameter  of  LOj. 

Then 


W  ~  uaP\\2L2(n)  <  kCq  Y  llu  ~  vi\\ and 


(13) 


\\V(u-uap)\\2L2{n)  <2 kY( 


ClWu-VjWl^ 

id,)2 


2.2.  Conditions  on  GFEM  data.  Our  main  results  will  involve  a  sequence  of 
GFEM-spaces  S„.  Let  {uij,  <pj,  be  a  single,  fixed  data  defining  a  GFEM- 

space  S,  as  in  the  previous  subsection,  and  let  E  :=  {ujj,(f>j,  Tj,  oj*},  where  w*  C  u>j. 
We  now  introduce  some  conditions  on  E  that  will  be  used  in  the  next  subsection 
to  formulate  our  assumptions  on  the  sequence  E„  =  {caf,  (ff,  'Ff,  defining 

GFEM-spaces  Sv,  v  G  Z+,  studied  in  this  paper. 
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Recall  that  co  is  star-shaped  with  respect  to  w*  C  co  if  for  every  x  £  co  and  every 
y  £  co* ,  the  segment  with  end  points  x  and  y  is  completely  contained  in  to. 

Condition  A(h).  We  have  that  =  U jLi^Oj  and  for  each  j  =  1,  2, . . . ,  N,  the  set 
LOj  is  open  of  diameter  dj  <  h  <  1  and  co*  C  Oj  is  an  open  ball  of  diameter  >  ah 
such  that  LOj  is  star-shaped  with  respect  to  to* . 

Condition  B.  The  family  {4>j}f—\  is  a  (k,  Cq,  C i,  . . . ,  Cm )  partition  of  unity. 

Condition  C.  For  each  j  =  1,  2, . . . ,  N,  the  space  \I >  j  contains  all  polynomials  of 
degree  A  and 

(14)  INI H‘(uj)  <  ^IMI 

for  any  w  £  'Sfj,  any  0  <  l  <  m,  and  any  ball  co*  C  u>j  of  diameter  >  ah. 

We  shall  need  to  use  “admissible”  subsets  of  fi,  a  class  of  subsets  that  we  now 
define. 

Definition  2.3.  Let  E  =  {tOj ,  <pj ,  4^,  co*}f=1  be  as  above  and  U  C  fl  be  an  open 
subset.  Denote  by  J(U)  the  set  of  those  indices  j  such  that  co*  C  U.  We  shall  say 
that  U  is  admissible  for  E  if,  for  all  j  =  1, . . . ,  N, 

(i)  c f>j  =  1  on  co*  and 

(ii)  J2jeJ(U)  $3  ~  1  011  U ■ 

We  can  now  formulate  our  last  condition. 

Condition  D.  The  domain  O  is  admissible  for  the  set  E  =  {u>j,  (j>j,  'bj,u)j}jL1 
defining  the  GFEM-space  S. 

Let  us  notice  that,  in  view  of  our  previous  conditions,  Condition  D  amounts  to 
the  fact  that  (j)j  =  1  on  cot  for  all  j  =  1, ... ,  N. 

2.3.  Assumptions.  We  are  now  ready  to  formulate  our  assumptions. 

(i)  We  assume  that  we  are  given  a  sequence  E„  =  {cof ,  (ff,  'LJ,  LOjl'}^=fv  v  £  Z+, 
defining  GFEM-spaces  Sv . 

(ii)  There  exists  constants  A ,  Cj,  a,  k,  A,  and  to  and  a  sequence  h v  — >  0, 
as  v  — >  oo,  such  that  E„  satisfies  Conditions  A(hv),  B,  C,  and  D  for  each 
v  €  Z_j^. 

The  constants  A ,  Cj,  a ,  k,  A,  and  to  will  be  called  structural  constants.  Note 
that  we  must  have  Nv  — >  oo  as  v  — »  oo.  For  simplicity,  we  shall  say  below  that 
the  sequence  Ej,  satisfies  the  Conditions  A-D  instead  of  saying  that  it  satisfies  the 
conditions  A (hv),  B,  C,  and  D. 

Let  us  recall  now  the  following  standard  lemma. 

Lemma  2.4.  Let  i pj  be  measurable  functions  defined  on  an  open  set  W.  Assume 
that  there  exists  an  integer  n  such  that  a  point  x  €  W  can  belong  to  no  more  than 
k  of  the  sets  supp (ipj).  Let  f  =  'Yhj‘4>j-  Then  there  exists  a  constant  C  >  0, 
depending  only  on  k,  such  that  ||/||^<(w)  <  W^jWh^w)- 

Proof.  The  inequality 

(15)  Ifli  +  d2  +  ■  •  ■  +  ojv; |2  <  M(|ai|“  +  | a>2 1 “  +  . . .  +  \om\2), 

with  k  <  M,  gives  the  desired  result.  □ 

We  can  now  establish  the  following  proposition. 
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Proposition  2.5.  Under  the  assumptions  of  this  subsection,  we  have  that  there 
exists  a  constant  B  >  0  such  that 

(16)  \w\h*{U)  <  -B/*v4|MU2(iO> 

for  all  0  <  t  <  m,  v  £  Z+,  w  £  Su,  and  U  C  £1  admissible  for  E^.  TTie  constant  B 
may  depend  only  on  the  structural  constants  A,  Cj,  k,  A,  cr,  and  m  (in  particular, 
it  is  independent  ofv,  u,  and  the  admissible  set  U). 

Proof.  Let  us  denote  h  =  hv,  for  simplicity.  Let  w  £  S.  Since  U  is  admissible,  we 
have  that  w  =  on  with  Wj  £  ’Lj,  the  sum  being  taken  over  all  j  such 

that  oj*  C  U .  Then  Lemma  2.4  and  Assumptions  A  -D  give 


\w\ HHU)  ^  CEE 

j  1=0 

<  Cj2J2cfh~2lB'2h 

j  1=0 


w 


^2K) 


<Ch-2t\\w\\lHU), 


where  all  the  structural  constants  C  above  depend  only  on  the  structural  constants 
A  Cjt  cr,  k,  A,  and  m.  This  proves  the  result.  □ 


We  obtained  right  away  that,  under  the  same  assumptions  as  those  in  the  above 
proposition,  that 

IMIff‘(C0  -  BKt\\w\\L2{u),  w  £  S. 

In  fact,  this  equation  is  equivalent  to  the  proposition.  Recall  that  the  constant  m 
is  the  fixed  integer  appearing  in  Assumptions  A-D. 


2.4.  Remarks.  We  now  include  a  few  simple  remarks  that  will  help  clarify  the 
above  conditions  and  assumptions. 

Remark  2.6.  Condition  A (h)  implies  that  the  diameters  of  uij  are  comparable  with 
h.  Indeed,  dj  :=  diamftCj j  >  ah. 


Remark  2.7.  The  explicit  inequalities  implied  by  Conditions  A  and  B  are 
(17)  II II we oo (Q)  <  Cih  l,  Z  =  0,l,...,m, 

by  Definition  2.1. 


Remark  2.8.  Condition  C  is  satisfied,  for  example,  if  4^  =  Q\,  the  set  of  polyno¬ 
mials  of  degree  at  most  A  >  1,  see  Section  6. 


Remark  2.9.  A  typical  example  of  an  admissible  set  is  obtained  as  follows.  Fix  a 
subset  J  of  indices  j  and  let  G  be  the  set  of  points  where  j  =  1-  Then  the 
interior  of  G  is  an  admissible  open  subset  of  fi. 


Remark  2.10.  If  U  C  O  is  an  admissible  open  set,  then  the  sets  u*  C  U,  must  be 
disjoint.  Moreover,  the  boundary  of  U  cannot  be  arbitrary,  see  Remark  6.7. 


10 


I  BABUSKA  AND  V.  NISTOR 


Remark  2.11.  Let,  for  each  fixed  j,  {wji}  be  a  basis  of  'F; .  Assume  that  Cl  is 
admissible.  Then  {ipjWji}  is  a  basis  of  the  GFEM-space  S.  This  is  an  important, 
non-trivial  consequence  that  is  not  always  satisfied,  see  for  example  [29,  30]. 

Finally,  let  IF  be  a  family  of  open  subsets  of  Cl.  We  shall  say  that  T  satisfies 
the  v-chain  condition  if  for  any  0  ^  Aq  <e  Au  C  O,  there  exist  Aq  <e  Bi  <g  B-2  <e 
. . .  <e  <g  Av,  with  all  Bj  £  T .  The  family  of  subsets  of  Cl  admissible  for 

satisfies  the  I'-cliain  condition  for  hv  small,  see  Lemma  6.5  (the  proof  of  this 
fact  requires  that  Cl  be  admissible  for  all  £„,  namely  Condition  D).  It  is  likely 
that  all  our  approximation  results  remain  true  if  one  replaces  Condition  D  with  the 
conditions  that  the  family  of  subsets  of  Cl  admissible  for  £„  satisfies  the  i'-chain 
condition  for  hv  small  enough. 

3.  Interior  estimates  for  the  GFEM 

From  now  on,  we  shall  assume  that  £„  =  {wj,  </>J, \I v  G  N,  is  a  sequence 
satisfying  all  assumptions  formulated  in  Subsection  2.3.  Also,  S„  will  be  the  result¬ 
ing  Generalized  Finite  Element  Space  of  Equation  (12)  associated  to  the  data  and 
hv  — >  0  will  be  the  corresponding  “small  parameters.”  In  this  section,  the  set  Cl 
will  be  a  bounded,  connected  open  subset  of  R"  (i.e.,  Cl  will  be  a  bounded  domain). 
We  shall  not  require  in  this  section  that  Cl  have  a  smooth  boundary. 

In  the  following,  we  shall  occasionally  drop  the  index  v,  in  order  not  to  overbur¬ 
den  the  notation.  For  instance,  we  shall  denote  S  =  Sv,  h  =  hv,  ujj  =  u>j,  and 
so  on,  when  the  index  v  is  understood.  Also,  recall  that  the  structural  constants  A, 
Cj,  a,  k,  A,  and  m  appearing  in  Conditions  A-D,  except  h  =  hv,  will  be  fixed  in 
what  follows. 

3.1.  ^-approximation.  We  shall  need  a  basic  result  on  the  approximation  of 
functions  in  Hs(Cl)  with  elements  in  the  GFEM-space  S,  extending  Theorem  2.2. 
Only  the  case  s  =  1  will  be  needed  in  this  section,  but  later  on  we  shall  also  need 
the  general  case. 

For  the  following  result,  we  shall  need  the  well  known  Bramble-Hilbert  lemma 
in  the  form  given  in  [9].  Let  us  recall  this  basic  result  in  the  form  that  we  need 
below,  as  well  as  the  relevant  definitions.  Let  to  C  be  a  bounded,  open  set  and 
let 

Pmax(w)  =  sup{  p ,  Cl  is  star-shaped  with  respect  to  a  ball  B  C  w  with  radius  p}. 

The  chunkiness  parameter  y(w)  of  uj  is  then  y(w)  :=  diam(w)//3max(w). 

If  /  is  a  smooth  function,  let 

QyJ,n(x)  =  f(y)+^2djf(y)(xj-yj)+...+  Y^  ~ — a!  =  or!  . . . aj, 

3  =  1  |a|=t 

be  the  Taylor  polynomial  of  /  at  y  of  degree  t.  If  B  C  u  is  an  open  ball,  then  Ql  f , 
the  Taylor  polynomial  of  degree  t  £  Z+  of  f  averaged  over  B  is  given  by 

(18)  Qtf{x)  =  f  Qv,f,n{x)<t>B{y)dy , 

Jb 

with  <f>B  £  Cf°{B)  a  function  with  integral  1.  (In  [9],  this  polynomial  is  called  the 
Taylor  polynomial  of  order  t  +  1  of  /  averaged  over  B.)  Note  that,  by  integration 
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by  parts,  we  can  extend  the  definition  of  Q*  f  to  /  €  L2(w).  We  assume  that  all 
the  functions  (j)B  are  affine  equivalent  to  a  fixed  given  function. 

We  shall  need  the  following  lemma  (Lemma  4.3.8  from  [9]). 

Theorem  3.1  (Bramble-Hilbert).  Let  u  be  an  open  set  with  chunkiness  parameter 
7 (to)  <  7.  Also,  let  B  C  u>  be  an  open  ball  with  radius  >  pmax(w)/2  such  that  w  is 
star-shaped  with  respect  to  B.  Let  Qt_1/  be  the  Taylor  polynomial  of  degree  t  —  1 
of  f  averaged  over  B.  Then 

1/  -  Qt~1f\H’(w)  <  C/,n,7diamHt_S|/lff‘(a;),  0  <  S  <  t. 

Here  Ct,nn  >  0  is  a  constant  depending  only  on  t,  n,  and  7. 

Recall  that  the  local  approximation  spaces  4^  contain  all  polynomials  of  degree 
<  A.  We  are  ready  now  to  prove  the  following  theorem. 

Theorem  3.2.  Let  U  C  O  be  an  admissible  subset  and  0<s<t<X  +  l,  s  <  m. 
Let  U'  =  U jU)j,  where  j  satisfies  u>*  C  U.  Then,  for  any  v  €  TT*  (L/') ,  there  exists 
w  €  S„  such  that 

||u  -  w\\Hs^u)  <  Cht~s  WvWh^u1) 

for  a  constant  C  that  depends  only  on  the  structural  constants,  and  is,  in  particular, 
independent  of  v,  v,  and  U. 

Let  us  notice  that,  by  taking  s  =  t  in  the  above  theorem,  we  immediately  obtain 
that,  using  the  same  notation,  that 

(!9)  IMIff‘(u)  <  CIMlHqu'). 

Also,  observe  that  O  =  U implies  that  U  =  U'  if  U  =  Cl. 

Proof.  Let  h  =  hv,  Wj  =  cuf,  and  so  on.  We  shall  use  the  notation  and  the  results 
from  [9] [Chapter  4]  introduced  before  Theorem  3.1.  Let  Wj(f)  =  Qt_1f  be  the 
Taylor  polynomial  of  degree  t  —  1  of  /  averaged  over  uj*,  for  u>*  C  U.  We  set 
Wj(f )  =  0  if  t  =  0  or  if  u>*  (fU.  Then,  by  the  Bramble-Hilbert  Lemma  recalled 
above  (Theorem  3.1),  we  have 

(20)  l/-^(/)I^K)<C'/it-s|/|Ht(a,.), 

with  a  constant  C  depending  only  on  s,  t,  and  a  (this  is  due  to  the  fact  that  ajj 
has  dimeter  <  hv  and  is  star-shaped  with  respect  to  the  ball  oj*  of  dimeter  >  ah, 
so  the  chunkiness  parameter  of  u>j  satisfies  7 (u>j)  <  a~l). 

Fix  now  v  €  iL4(Z7')  arbitrary.  Also,  let  Wj  =  Wj(v )  G  4^  and  w  =  G  S, 

the  sum  being  taken  over  all  indices  such  that  u>*  C  U.  Then,  using  also  Condition 
B,  Lemma  2.4,  and  y~b  (f>j  =  1  on  U ,  we  obtain 


(21) 


I  H-(U) 


<  4>j(y  - 


j  i—1 


wj)\Hr-i(  Wj 


<cEEc*rt2t_2r+2i<*(^) 


j  i—l 


<  CKh2t~2r\V \2Ht(ul) 


Summing  over  0  <  r  <  s  and  using  0  <  h  <  1  gives  the  desired  result. 


□ 
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Let  us  also  record,  for  further  use,  the  following  well  known  Poincare— Friedrichs 
inequality  [9],  Lemma  (4.3.8).  (See  also  [9],  Lemma  (4.3.14),  and  [10],  Equation 
(2.2),  Theorem  14.1,  and  Theorem  15.3.,  or  [11,  31].)  The  precise  statement  that 
we  need  is  the  following. 

Theorem  3.3.  Assume  to  is  an  open  set  of  diameter  diam(w)  star-shaped  with 
respect  to  a  ball  u>*  C  u>  of  diameter  >  udiam(w).  Let  0  <  if  <  1  be  a  measurable 
function  on  u>,  if  =  1  on  u>* ,  and  v  =  f  v(x)ip{x)dx  /  J ip(x)dx  be  the  weighted 
average  of  v  overuj.  Then  there  exists  a  constant  Cp  that  depends  only  on  a  and 
the  dimension  n,  but  not  on  ip  or  u>,  such  that 

(22)  \\v-v\\l2(u)  <  Cpdiam(w)Mffi(w), 

for  all  v  €  iL1(w). 

When  ip  =  4>b  from  Equation  (18),  the  result  above  reduces  to  the  Bramble- 
Hilbert  lemma.  Except  for  the  fact  that  Cp  depends  only  on  n  and  h,  this  result  is 
well  known  when  v  is  the  average  of  v  over  to.  The  more  general  form  of  the  above 
theorem  may  be  useful  when  checking  that  the  assumptions  of  Subsection  2.3  are 
satisfied  for  a  suitable  sequence  Sv  of  GFEM-spaces. 

Throughout  this  paper,  C  will  denote  a  generic  constant  that  depends  only  on 
the  dimension  n. 


Proof.  Using  a  dilation  and  the  homogeneity  properties  of  the  norms  in  the  state¬ 
ment,  we  see  that  we  can  reduce  to  the  case  h  =  <j~1  and  o>*  =  B i,  the  ball  of 
radius  1  centered  at  the  origin.  Then  u>  C  B2h  =  B2(7-x.  Let  (p  be  the  function 
appearing  in  the  definition  of  the  averaged  Taylor  polynomial  Qk  f,  Equation  (18). 
Thus  (p  has  support  in  the  unit  ball  u>*  =  B\ . 

Recall  that  C  denotes  a  generic  constant  that  depends  only  on  the  dimension  n. 
Let  v  £  H 1(w)  be  arbitrary.  We  subtract  from  v  various  constants  to  obtain,  using 
also  the  Bramble-Hilbert  lemmma,  Theorem  3.1,  for  l  =  1  and  k  =  1, 


\\v  -  V  ||l2(w)  <  ||u  -  /  (p(x)v{x)dx\\L2^  +  \\  /  cp(x)v{x)dx  ~  v\\L2^ 

J  U)  J  UJ 

<  C,|M|ff1(w)  +  ll  /  (p{x)v(x)dx-v\\L2(U})  <  C(\\v\\  + 1  /  <P{x)v(x)dx-v\), 

J  UJ  J  UJ 


where  by  Jw  cp(x)v{x)dx  —  v  we  mean  the  constant  function  on  oj  with  this  value 
(so  its  L2  norm  is  a  multiple  of  this  constant,  and  in  our  case  this  multiple  can  be 
bounded  by  a  constant  depending  only  on  the  dimension  n ). 

Let  a  :=  J  ip(x)dx  and  g  =  <p  —  a~lip,  so  that,  in  particular,  f^g(x)dx  =  0. 
Since  v  =  a^1  J '  v(x)ip(x)dx,  it  is  enough  to  show  that 


(23) 


g(x)v(x)dx  <  C\v\h^(, 


Let  i S™-1  =  dBi  denote  the  unit  sphere  of  radius  1  in  R"  (the  boundary  of  the 
unit  ball  B i).  Let  cp\  have  support  in  the  closure  of  Bi  =  u>* ,  be  constant  on  the  ray 
Rx7  fl  Bi,  and  satisfy  /0°° {(pi{rx')  —  a~1ip(rx'))rn~1dr  =  0,  almost  everywhere  in 
x1  £  S'"  .  (Note  that  the  function  to  be  integrated  vanishes  if  r  >  2<r”1.)  Then  <pi 

is  bounded  by  a  constant  that  depends  only  on  the  dimension.  Let  g\  =  <pi  —  a~lrip 
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and  let  G{tx')  =  /0*  g\(rx')rn  1dr ,  so  that  \G(tx')\  <  Ctn.  An  integration  by  parts 
then  shows  that 


gi  ( x)v{x)dx 


G(x)drv  (cc)  |a;| 1  ndx, 


J  UJ  J  UJ 

where  dr  is  the  derivative  in  the  radial  direction.  Since  G(x)\x\1~n  is  bounded  by 
a  constant  that  depends  only  on  the  dimension  n ,  we  obtain  that 


gi(x)v(x)dx  <C 


which  is  the  desired  Equation  (23),  but  with  g  replaced  by  g\.  It  is  therefore  enough 
to  prove 


[gi(x)  -  g{x)]v{x)dx 


gi(x)  -  g{x)\v{x)dx  <C\v\Hi(u) 


where  the  first  equality  is  due  to  the  fact  that  g\  —  g  =  </>i  —  f>  has  support  in  B\ . 

Let  du  denote  the  derivative  in  the  direction  of  the  outer  normal  to  80.  Since 
fB  [g-\  (x)  —  g(x)]dx  =  0,  we  can  find  V  be  such  that  AV  =  g\  —  g  and  dvV  =  0 

on  the  boundary  of  B\.  Moreover,  ||  W||£2  <  C\\g\  —  g\\i,-2  <  C,  where,  we  recall, 
C  is  a  generic  constant  that  may  depend  only  on  the  dimension  n.  Then 


'  [9i(x)  -  g(x)\v(x)dx 

UJ 

The  proof  is  now  complete. 


VE  •  'S/v(x)dx  <  C\v\h^(u,)- 


□ 


For  k  =  1,  we  shall  need  the  following  consequence  of  Theorem  3.2,  which  re¬ 
places  Assumption  9.5  of  [33]  and  does  not  require  the  open  sets  involved,  except 
O,  to  be  admissible.  Define 

(24)  S<(n):=SvnCc(Q). 

That  is,  S<(Q)  consists  of  the  elements  of  the  GFEM-space  S„  with  compact 
support  inside  fl. 

Recall  that  A  <g  B  means  that  the  closure  of  A  is  a  compact  set  contained  in 
the  interior  of  B  (i.e.,  A  is  a  relatively  compact  subset  of  B).  Also,  recall  that  the 
family  T  of  subsets  of  admissible  for  E„  satisfies  the  i/-chain  condition  for  hv 
small  enough,  see  Subsection  2.4  and  Lemma  6.5. 

Our  main  goal  in  this  section  is  to  prove  Theorem  3.12. 


Proposition  3.4.  Let  U  <e  Oi  C  LI  subsets  of  O  and  9  be  the  distance  from  8U 
to  dfli.  Assume  that  U  is  admissible.  Then  there  exists  C  >  0,  independent  of  9, 
v,  and  U ,  with  the  following  property.  For  any  u  €  H2(Ll)  with  support  in  U,  there 
exists  w  £  such  that 


if  K  <  0- 


||u  —  HI-h-i^)  —  C,^llull/r2(ni)> 


Proof.  Choose  Wj  and  w  as  in  the  proof  of  Theorem  3.2.  We  shall  continue  to 
use  the  notation  of  that  Theorem.  In  particular,  h  =  hv.  If  h  <  9,  then  Wj  =  0 
unless  ajj  intersects  U,  which  gives  that  the  closure  of  uij  is  completely  contained 
in  Hi.  In  particular,  U'  C  ff|  and  the  support  of  w  constructed  above  is  compact 
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and  contained  in  Q) .  Then  we  can  replace  U  with  U'  in  Equation  (21)  and  we  thus 
obtain 

||u  —  w||_f/i(Oi)  =  ||w  —  w\\ hhu')  <  Ch\\u\\H^(u')  =  C'/i||u||/J-2(q1). 

This  completes  the  proof.  □ 

Remark  3.5.  By  taking  (7=1  and  w  =  0  for  h>  9,  we  obtain 
II w  —  <  C9  1/iI/||u||ff2(Q1), 

for  all  v  (not  just  for  hv  <  9). 

3.2.  The  super-approximation  property.  The  assumptions  of  Subsection  2.3 
continue  to  remain  valid,  in  particular,  all  the  structural  constants  below  will  be 
independent  on  v.  Also,  recall  that  we  have  considered  in  the  Introduction  the 
bilinear  form 

(25)  B(u,v)  :=  I  Vu-Vvdx,  u,v£H1( £2). 

Jn 

Our  approach  follows  the  approach  from  [21],  as  presented  in  [33]  [Section  9].  See 
[8,  20,  23,  24]  for  related  results  on  approximation  in  the  “sup”-norm. 

Lemma  3.6.  Let  p  be  a  smooth  function  on  u>j  and  id  £  fj.  Then  there  exists 
w  £  ^ j  such  that 

(26)  || pw  -  w\\m{u>3)  <  Ch\\p\\W2,ao^.)\\w\\Hi{uij)% 

where,  we  recall,  <7  >  0  may  depend  only  on  the  dimension  n  (in  particular,  it  is 
independent  ofw,  p,  andj). 

Proof.  We  shall  use  the  inner  product  induced  from  Hl{uij).  Let  p  £ 
be  given. 

To  prove  the  lemma,  we  shall  assume  first  that  w  £  Tj  is  a  constant.  Let  L 
be  the  degree  one  Taylor  polynomial  approximation  of  p  at  the  center  of  the  ball 
to*.  Then  L  €  fj,  because  first  order  polynomials  are  in  Tj  (Condition  C).  Since 
h  <  1,  we  obtain 

\\p  -  L\\wi <  Ch\\p\\w^.y 

Choose  w  =  Lw.  Then 

\\pw  -  w\\Hi(Uj)  =  \\pw  -  Lw\\hi{uj)  <  ||p  -  L\\wi,cc{uj)  HujIIhi^) 

<  Ch\\p\\W2,°o(u.)\\w\\Hi(u.y 

Assume  now  that  w  £  is  such  that  (w,  1)  =  0,  that  is,  w  is  orthogonal  in 
i71(wj)  to  the  subspace  generated  by  constants.  We  then  write 

P= P+P , 

where  p  is  a  constant  function  (say  the  value  of  p  at  the  center  of  u>*)  and 

(27)  noo(uj)<dh\\Vp\\Loo(uj). 

We  shall  choose  then  w  =  pw  £  Tj,  which  makes  sense  since  Tj  is  a  vector  space. 
Then 

\\pw  -  w\\Hiiujj)  =  \\p*w\\Hi{u}j)  <  <7||Vp*||ioo(w.)||u;||i2(a,.) 

+  C||P*l|L~(a;dllU'll«1(wd  -  C'^ll^7PlU00(w3)llU,Hffl(w3)’ 
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where  in  the  last  step  we  have  used  the  Poincare-Friedrichs  inequality  for  uij  (The¬ 
orem  3.3)  to  estimate  ||w||l2(w.)  and  Equation  (27)  above  to  estimate  the  second 
term. 

For  a  general  w  £  \E '  j,  we  decompose  w  =  w±  +  w2  with  W\  a  constant  and  w2 
orthogonal  to  the  space  of  constants  and  choose  w-\  and  w2  as  above.  Then 

II pw  -wi-  <  \\pwi  -  wiWhi^j)  +  \\pw2  - 

<  Ch\\p\\W2,oo(u,.')\\w1\\Hi(u,.')  +  IU2'0°(uJj)  ll^72 1|  H1  {ujj) 

<  Ch\\p\\wz.°° (u^WwWmiwj)- 

The  lemma  is  now  proved.  □ 

Recall  that  we  denote  A  <g  B  if  A,  the  closure  of  A  in  R2,  is  a  compact  subset 
of  the  interior  of  B.  Also,  recall  from  Equation  (24)  that  Sf'(A')  denotes  the  set  of 
elements  in  Sj,  with  compact  support  contained  in  A,  for  any  open  subset  A  C  Cl. 

An  important  technical  step  in  our  proof  of  the  Theorem  3.12  is  the  following 
“super-approximation”  result. 

Proposition  3.7.  Let  p  £  W2’°°( LI)  and  iu  £  Sv.  Then  there  exists  w  £  Sv  such 
that 

(28)  \\pw  —  ■u’Hffi(n)  <  C,/i^||p||w2>~(n)||w||_f/i(n), 

where  C  is  independent  of  v.  If  w  has  support,  in  fii  <s  Cl2  and  0  is  the  distance 
from  dLl\  to  dd2,  then  w  £  S<(fl2)  for  h„  <  9. 

As  explained  above,  the  constant  C  may  depend  on  the  structural  constants  A, 
Cj,  k,  to,  <7,  and  A,  but  is  independent  of  h  =  hv  and  of  the  number  Nv  of  sets 
{coj}.  In  particular,  it  is  independent  of  the  GFEM-space  Sv  used. 

Proof.  Let  h  =  hv  and 

N 

(29)  w  =  < i>jWj  £  Sv,  Wj  £  Tj. 

i-J 

Let  Wj  be  the  orthogonal  projection  of  pwj  onto  'Fj  in  the  inner  product  of  H1[wj'). 
Lemma  3.6  then  shows  that 

(30)  || pwj  -  Wj\\H i(Wj.)  <  Ch\\p\\W2,°o(Q,)\\wj\\Hi(W:j). 

Moreover,  we  have  that  f  ( pwj  —  Wj)dx  =  0  because  the  constant  functions  are  in 
Tj  and  pwj  —  uij  is  orthogonal  to  'Pj. 

Let  w  :=  Then  HV^Hloo^)  <  C\/h  by  Equation  (11)  and 

II pwj  -WjWmuj)  <  CPh\\pWj  <  CpCh2\\p\\W2,^^\\wj\\Hi(u.), 

by  the  Poincare-Friedrichs  inequality  (Theorem  3.3),  and  hence 
N  N 

\\pw-w\\2HHn)  =  \\'52<t>j(pWj-Wj)\\HHn)  <  CY1 

3=1  3=1 

N 

+  l|V^|||ooK)||p^  -  Wj\\2L 2(Wj))  <  C'/i2||p||^2>0o(Q)  y]  \\wj\\2Hi(uj)i 

3=1 
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where  for  the  first  inequality  we  have  used  also  Lemma  2.4.  The  result  will  follow 
now  if  we  can  prove  that  ^  CIMIffqnp  for  any  w  =  EyLi  fawj, 

as  above  and  C  a  constant  independent  of  5  =  S„.  Indeed,  since  (fij  =  1  on  ui*.  by 
Condition  D,  we  have  w  =  Wj  on  lu*  ,  and  hence 

N  N 

IMlirqfi)  ^  \\wj\\m(u>])  ^  ^  2  Ill'll 

J'= 1  3  = 1 


by  Condition  C  (A  is  the  constant  appearing  in  that  assumption). 

The  proof  of  the  last  part  is  completed  as  in  Proposition  3.4.  □ 


3.3.  Estimates  on  “discrete— harmonic”  functions.  We  shall  also  need  the 
following  “inverse  property,”  which  is  somewhat  similar  to  Assumption  A. 3.  in  [21] 
or  Assumption  9.2  in  [33]. 

The  rest  of  this  section  follows  closely  the  approach  in  the  paper  of  Nitsche 
and  Schatz  [21],  relying  also  from  the  survey  paper  [33]  (which  in  turn  is  based 
on  the  paper  by  Nitsche  and  Schatz).  There  are,  however,  some  differences  in  the 
assumptions  that  we  are  using,  so  we  include  complete  proofs  for  the  convenience 
of  the  reader.  For  instance,  the  following  lemma,  Lemma  3.8,  plays  the  role  of 
Assumption  A. 3.  in  the  Nitsche-Schatz  article  [21],  respectively,  of  the  Assumption 
9.2  (Inverse  assumption)  in  Wahlbin’s  article.  Also,  the  following  lemma  is  an 
analog  of  Lemma  5.2  of  [21],  respectively,  of  Lemma  9.1  of  [33]. 

Lemma  3.8.  There  exists  a  constant  C  >  0,  depending  only  on  the  structural 
constants,  such  that 

IM|l2(e/)  <  Ch-j\\w\\H-i(U), 

for  all  0  <  j  <  m,  is  £  Z+,  w  £  Sv,  and  U  C  O  admissible  for 


Proof.  Let  h  =  hv.  We  have 


IM|i?-j(fi)=  sup 
0 


\(w^)\ 


> 


by  Proposition  2.5. 


(w,  w) 

IMI-FD(n) 


IMI/tf  (n) 


>C  1hJ\\w\\L2^n), 


□ 


We  now  prove  the  following  crucial  lemma. 

Lemma  3.9.  Let  A  <s  A'  <e  O  be  open  subsets.  Then  there  exists  C  >  0  with  the 
following  property.  If  w  £  Su  and 

(31)  B(w,x )  =  0,  for  all  \  G  S<(A'), 

then,  for  h„  small  enough,  ||u>||.ffi(A)  <  C||wI|l2(A')j  C  depending  on  the  dis¬ 
tance  6  from  d A  and  d A! ,  but.  not  on  is  or  w. 

Proof.  The  proof  is  very  similar  to  the  one  in  [21,  33],  using  Proposition  2.5  in 
place  of  Condition  A. 3,  respectively  Assumption  9.2  (“Inverse  assumption”),  and 
Proposition  3.7  in  place  of  Assumption  A. 2,  respectively  Assumption  9.1  (“Super¬ 
approximation”),  of  [21],  respectively  [33]. 

Let  h  =  hv  be  small  enough.  Let  us  chose  open  sets  A  <e  A0  <e  A\  <g  A!  with  A\ 
admissible  for  Ej,  and  with  the  distances  between  the  boundaries  comparable  with 
9.  Also,  let  uj  £  Cf°(A0),  with  to  =  1  on  A  and  w  >  0  on  A0  such  that  all  the  norms 
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||w||vufc,oo  are  bounded  by  a  constant  depending  only  on  9.  Then,  by  Equation  (31), 
we  obtain 

II^HIl2(A)  —  (VuywVre)  =  (Vw,V  (ojw))  —  (Vw,V(w)w) 

=  {S7w,S7{low  —  ip))  +  ^(ic,  (Au)u>), 

where  the  inner  products  are  in  L2(Ao)  and  ip  £  S)f(Ao).  Proposition  3.7  then 
gives 

ll^wlli2(A)  <  C^IMIw2.~(fi)|M|lP(Ao)  +  ^1Mll2(A0)> 

which,  in  turn,  implies 

(32)  IMIffi(A)  <  C(h1/2\\w\\Hi(Ao)  +  |MU2(A0)), 

where  all  the  constants  depend  only  on  9.  We  now  repeat  the  argument  for  A$  <g 
A-\  (e  Q  (and  A  replaced  by  A0  and  A0  replaced  by  Af),  which  gives 

(33)  IMIffTAo)  <  C(h1/2\\w\\HHAl)  +  ||u;||L2(j4l)). 

Combining  Equations  (32)  and  (33)  and  using  also  h  <  1,  we  obtain 

(34)  Ikllffi(A)  <  C(h\\w\\Hi{Al)  +  C\\w\\LHAl)). 

Then  we  use  Proposition  2.5  to  obtain  that  h\\w\\jj1(A1)  <  C||u;IIl2(Ai)i  with  C 
independent  of  w  and  v  £  Z+.  This  finally  gives 

(35)  lkl|j?i(A)  <  C'I|w||l2(a1)  <  C\\w\\L2(A'), 

with  a  constant  depending  only  on  9.  The  proof  is  now  complete.  □ 

We  shall  need  the  following  simple  estimate. 

Lemma  3.10.  Let  $(;r)  =  log  \x\  if  n  —  2,  $(a:)  =  \x\2~n ,  if  n  ^  2.  Let  U  cLl  be 
an  open  subset.  Then  there  exists  C  >  0,  independent  of  U,  such  that 

$  *  u(x)  :=  /  $(x  —  y)u(y)dy 

Ju 

satisfies 

||$  *  u\\Hi+2(jj)  <  Cl||u||H((!7), 

for  any  l  £  Z  and  any  u  £  Cf°{U).  The  constant  C  will  depend  on  fi,  however. 

Proof.  Fix  U  C  Ll  and  let  u  £  Cf°(U)  C  Cf°(fl)  and  v  =  $  *  u.  We  have  Av  =  cnu 
for  some  constant  cn  depending  only  on  the  dimension  n  (this  well  known  fact  is 
proved,  for  example,  in  [11]).  The  generic  constants  below,  denoted  C,  are  allowed 
to  depend  only  on  R  and  the  dimension  n. 

Let  O  =  Br(0)  be  the  ball  of  radius  R  >  3diam(fl)  centered  at  the  origin,  where 
diam(fl)  is  the  diameter  of  O.  We  shall  assume  that  R  is  very  large.  In  particular, 
we  shall  assume  that  Q  C  BR/3( 0).  Let  rj  :  [0,oo)  -»  1  be  a  smooth  function  such 
that  rj (t)  =  1  if  t  <  diam(fi)  and  r?(f)  =  0  if  t  >  2diam(fl).  Let  $i(cr)  =  77(|x|)$(x) 
and  define  v  =  $i  *  u. 

Then  v(x)  vanishes  if  the  distance  from  x  to  f 1  is  greater  than  2  diam(fl)  <  2 R, 
so  v  =  0  on  dO  =  dBR{ 0).  Also,  $(x)  =  $i(x)  for  |x|  <  diam(fi),  and  hence 
v(x)  =  $  *  u(x)  for  x  £  O.  Moreover,  Av  =  (A$i)  *  u  =  cnu  +  <p  *  u,  where 
p(x)  =  A$i(x)  for  x  ^  0  and  <p(x)  =  0  if  |x|  <  diam(fl).  Since  p  £  L1(M”), 
we  have  that  \\p  *  and  hence  UAv^^n)  < 

Finally,  since  v  =  0  on  dO ,  we  obtain  that  IMI/fi+zm)  <  C'||tt||ffqRr*),  by  standard 


18 


I  BABUSKA  AND  V.  NISTOR 


estimates  on  elliptic  boundary  value  problems  (see  [11]  or  [31],  for  instance).  This 
proves  that 

11$  *  u\u\\h‘+2{U)  =  IMt/|l//!+2(E/)  <  CIMIffipt")  =  C\\u\\h1(U)- 
The  proof  is  complete.  □ 

An  alternative,  shorter  proof  of  this  lemma  can  be  obtained  using  the  fact  that, 
for  any  smooth,  compactly  supported  function,  u  — >  w(d>*u)  is  a  pseudodifferential 
operator  of  order  —2  with  compact  distribution  kernel  [32]. 

We  define 

(36)  IMIg-s([/)  =  SUP  iifiT’  ^  <  lM|g-(t/),  Cf°  (U) 

for  any  open  set  U,  any  u  £  L2{U),  and  any  seN.  We  define  HqS(U)  to  be  the 
completion  of  L2(U )  in  the  norm  IMI#- »(m-  Then  HqS(U),  s  €  N,  identifies  with 
the  dual  of  Hq(U). 

Let  9  be  the  distance  from  d A  and  dA' ,  as  before.  We  now  prove  the  following 
lemma. 

Lemma  3.11.  We  keep  the  notation  and  assumption  of  Lemma  3.9.  In  particular, 
we  assume  that  w  £  Sv  satisfies  Equation  (31).  Then,  for  h „  small  enough, 

(37)  IMIl2(A)  <  C'lMl.w— m(A')> 

where  C  is  a  constant  depending  9  and  the  structural  constants,  but  not  on  v. 
Combining  Lemmas  3.9  and  3.11,  we  obtain 

(38)  IMIffbA)  <  C\\w\\H-n,(A'), 

for  /i„  small  enough  and  any  w  £  Sv  satisfying  the  assumptions  of  Lemma  3.9. 

Proof.  Let  h  =  hu  and  A  <g  Bq  <s  B\  <e  A'  be  open  sets  with  Bq  admissible.  Note 
that  Lemma  3.9  gives 

(39)  IMIffT-Bi)  <  C\\w\\L2(AQ. 

For  any  v  £  Cf°(A),  let  V  :=  cn&  *  v  £  Hl+2(B i),  where  cn  is  chosen  such  that 
AC  =  v  (see  [11]).  Lemma  3.10  then  gives 

(40)  W\\h1+2{b1)  <  C||u||ff!(Bl)  =  C\\v\\Hi(A),  l  £  Z+, 

for  some  constant  C  that  depends  only  on  H. 

Let  u>  £  Cf°(B0)  with  w  =  lon  A.  Since  luV  £  Cf°(B0)  and  Bq  is  admissible,  we 
know  from  Proposition  3.4  that,  for  h  small  enough,  there  exists  \  €  Sf  (B\)  such 
that 

(41)  || uV  —  xI|hi(B!)  <  Ch\\u)V\\H2(B1)  <  Ch\\V\\H^(B1)  <  Ch\\v\\L2^A)- 

Then,  still  assuming  that  v  £  C£°(A)  is  arbitrary  and  taking  into  account  also 
Equation  (31),  we  obtain 


(w,v)a  =  (u>w,v)a  =  /  u>wAVdx  =  /  uiwAVdx  =  —  /  V(wtu)  ■  Wdx 

J  A  J  Bq  J  Bq 

=  —  w(2\7u)  ■  X/V+V Aufjdx—  /  Vie  ■  V(cjV  —  x)dx, 

J Bq  J  B\ 
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for  any  y  G  S< (Bi)  C  S^(A'),  where  the  inner  products  are  calculated  on  the 
indicated  sets.  Then,  by  combining  Equations  (39),  (40),  and  (41),  as  well  as 
Lemma  3.9,  we  obtain,  for  all  l  >  0, 

|(w,^)a|  <  C|H|ff-<-i(Bo)||E||_fP+2(Bl)  +  C/i||w||Bi(Bl)||V'||B2(Bl) 

<  C'dkllff-'-^A')  +  h\\w\\L^A'))\\v\\H‘(A), 

and  hence 

(42)  lk||H-i(A)  <  C(\\w\\H-i-i{AI)  +  h\\w\\L2{AI)). 

Next,  let  us  choose  a  sequence  of  open  sets  A  <g  Bi  <g  B2  <e  . . .  <s  Bm  <g  A'  with 
the  distances  between  the  boundaries  comparable  with  9.  Changing  notation  and 
iterating  Equation  (42),  using  also  h  <  1,  we  obtain 

(43)  |MU2(A)  <  C(\\w\\H-i{Bi)  +  h\\w\\L2{Bl)) 

<  C(\\w\\h-2(B2)  +  h\\wh 2(B2))  < 

<  C(\\W\\H-™(Bm)  +  h\\W\ U2(B™))- 

We  now  repeat  the  above  reasoning.  We  change  notation  again,  so  that,  this 
time,  Brn  becomes  B\,  then  we  chose  as  before  a  sequence  of  open  sets 

(44)  A  <g  Bx  <e  B2  <e  . . .  <e  Bm  <g  A'. 

We  can  assume  that  the  distances  between  the  boundaries  are  comparable  with 
6.  Also,  we  can  assume  that  Bm  is  admissible,  since  the  family  of  admissible  sets 
satisfies  the  //-chain  condition  for  h  small  enough,  see  Lemma  6.5.  Then  we  iterate 
Equation  (43)  using  again  h  <  1,  and  obtain, 

IMU2(A)  <  C(\\w\\H-m{Bl)  +  h\\w\\L2{Bl)) 

<  C'dkll  H-™(B2)  +  ^2|IwIU2(B2))  <  •  •  •  <  C'dkll +  hm\\w\ \L*(Bm)) 

<  C^IMI ^  C\\w\\H-m(Bm), 

where  at  the  end  we  have  used  the  inverse  property  ||u’||l2(j/)  < 

for  any  admissible  open  subset  U  C  O,  (see  Lemma  3.8).  The  proof  is  complete.  □ 

3.4.  The  interior  error  estimate.  The  following  result,  the  main  result  of  this 
section,  is  an  analog  of  [21]  [Theorem  5.1]  and  of  [33]  [Theorem  9.2]. 

Let  Sv  be  the  GFEM-spaces  associated  to  the  data  uf,  4't',  u/*1'} 

satisfying  the  Assumptions  of  Subsectio  2.3. 

Theorem  3.12.  Let  A  <s  B  C  fi  be  open  subsets.  Then  there  exists  C  >  0 
with  the  following  property.  If  u  G  H1(  12)  and  u „  G  Sv,  v  G  Z+,  are  such  that 
B(u  —  uu,x)  =  0  Sor  X  G  ^(fi),  then  for  hv  small  enough, 

||«-  uv\\m(A)  <  inf  ||u-xl|ffi(B)  +  \W  ~  uv\\H-m(B^j . 

The  constant  C  depends  only  on  the  distance  from  dA  to  dB  and  the  structural 
constants,  but  not  on  v  G  Z+. 

For  example,  uu  in  the  above  theorem  could  be  the  GFEM-approximation  of  u 
(see  Equation  (54)),  but  our  assumptions  on  uv  are  in  fact  somewhat  weaker. 
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Proof.  Let  A  <g  A\  <g  A2  <s  B  c  fi.  Choose  u>  G  C£°(A2)  such  that  u  =  1  on  A\. 
Let  Pi  be  the  -fP(fl)  orthogonal  projection  onto  Sf  (A-| )  C  S  C  iP(fi).  Then  on 

(45)  u  —  uv  =  (uu  —  Pi(iou))  +  (Pi(unt)  —  it„). 

Then,  by  the  general  properties  of  orthogonal  projections,  we  have 

(46)  \\uju  —  Pi(wu)||j?i(Q)  <  ||wM||ffi(Q)  <  C\\u\\Hi(By 
Hence 

(47)  ||wm  -  Pi(oJu)\\H~m^Al)  <  \\oju  -  Pi(ww)||ffi(fi)  <  C'HuHhi(b). 

Let  ic  =  Pi(um)  —  uv.  Then  B(w,  x)  =  B(u>u  —  u„,  \)  =  B(u  —  uu,  x)  =  0,  for  all 
X  G  iSjf  (Ai),  and  hence  w  satisfies  the  assumptions  of  Lemmas  3.9  and  3.11.  From 
this,  using  also  Equations  (38),  (45),  and  (47),  we  obtain 

(48)  IHIffi(A)  <  C\\w\\h-™(Ai)  <  II WM  -  Pl(wu)||H-m(Al)  +  ||u  -  uv\\H-m{M) 

<  \\u\\hi(B)  +  ||«  -  Uv\\H~m(Al ). 

Equations  (45-48)  then  give 

\\u- uv\\Hi{A)  <  \\um-  Pi(uju)\\hi{a)  +  ||w||bi(a)  by  (45) 

<  C\\u\\Hi(B)  +  ||w  -  u„||ff-m(Al)  by  (46)-(48). 

The  desired  result  follows  by  replacing  u  and  uu  with  u  —  x  and,  respectively,  uu  —  x, 
with  %  in  S  =  Sv.  □ 


4.  Discrete  solutions 

We  assume  in  this  and  next  section  that  Q  =  O,  a  smooth,  bounded  domain. 
We  discuss  in  this  section  a  variational  formulation  of  the  boundary  value  problem 
(1).  We  also  prove  the  existence,  uniqueness  of  the  solution  of  this  problem  as  well 
as  the  estimate  (2).  (So  we  establish  the  well  posedness  of  this  boundary  value 
problem.)  The  discrete  solutions,  or  GFEM-approximations  to  u  are  also  defined 
using  the  weak  formulation  of  (1).  We  begin  with  the  definition  of  the  discrete 
solutions.  The  results  extend  to  the  case  of  g  with  low  regularity  the  usual  results 
and  definitions  of  the  weak  solution  and  its  discretization,  thus,  the  reader  willing 
to  accept  that  the  theory  is  essentially  the  same  as  in  the  case  g  G  L2{dO),  can 
skip  this  section  at  a  first  reading. 

Let  us  fix,  in  this  section,  a  finite  dimensional  subspace  S  C  Hm(0)  containing 
all  constant  functions.  In  applications,  S  =  Sv,  for  some  v  G  Z+,  but  that  is  not 
required  in  this  section. 

4.1.  Discrete  solutions.  Consider  the  usual  bilinear  form 
(49)  P(0,V>):=  /  V0  •  V0  dx, 

Jo 

where  <j),ip  G  Hl(<D).  Assume,  for  the  purpose  of  this  discussion,  that  H  =  O  is 
smooth,  so  that  the  boundary  value  problem  (1)  is  defined. 

Recall  that  we  write  (it,  v)qo  '■=  u(v)  G  C  for  any  u  G  P1_fc(0)  and  v  G  iLfc_1(0) 
for  the  “value  of  u  evaluated  at  v.” 
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Lemma  4.1.  Assume  ( g ,  l)go  =  0  and  g  £  H  1/2  k(dO),  k  <  m  —  1.  Then  there 
exists  a  unique  u$  G  S  such  that  ( us ,  1)  =  0  and 

(50)  B(us,vs)  =  {g,vs\do)dO,  for  all  vs  £  S, 

for  all  vs  £  S.  (Recall  that  S  C  Hm(0).) 

Note  that  S  C  Hm(0)  and  g  £  H~1J2~k(dO),  k  <  m  —  1,  guarantees  that 
{g,vs\do)dO  is  defined. 

Proof.  Let  So  be  the  subspace  of  the  GFEM-space  S  consisting  of  functions  Xo  £  S 
with  (xo>  1)  =  Jq  Xo(x)dx  =  0.  The  bilinear  form  B  is  non-degenerate  on  So  (that 
is,  if  (p  €  So  is  such  that  B(<p,ip)  =  0  for  all  ip  £  So,  then  <p  =  0).  This  gives  by 
standard  linear  algebra  the  existence  of  a  unique  us  £  So  such  that  Equation  (50) 
is  satisfied  for  all  vs  £  S0.  □ 

Definition  4.2.  If  us  £  S  is  as  in  the  above  Lemma,  then  we  shall  say  that  us  is 
the  GFEM-approximation  of  the  solution  of  Equation  (1). 

Similarly,  if  g  £  id1/2(0)  and  (g,  l)oo  =  Jg0  dix)dS(x)  =  0,  then  the  solution  u 
of  Equation  (1)  is  in  H2((D),  ( u ,  1)  =  0  and, 

(51)  B(u,v)  =  (g,v\d0)do, 

for  all  v  smooth  enough.  Moreover,  u  is  uniquely  determined  by  these  conditions. 
This  is  of  course,  nothing  but  the  weak  formulation  of  the  boundary  value  problem 
of  Equation  (1).  We  now  extend  the  formulation  to  the  case  when  the  data  g  is  a 
distribution. 

Lemma  4.1  states  only  the  existence  and  uniqueness  of  the  discrete  solution  us 
of  Equation  (1).  It  claims  nothing  about  the  relation  between  us  and  the  exact 
solution  u  of  Equation  (1).  This  will  be  discussed  in  the  remaining  of  this  section 
and  in  the  following  section. 

4.2.  The  weak  solution.  Let  us  begin  with  a  remark  that  will  justify  the  following 
constructions. 

Remark  4.3.  When  trying  to  extend  the  weak  formulation  of  Equation  (1)  to  g 
with  lower  regularity  ( i.e g  a  distribution),  we  face  the  following  difficulty.  Let 
v  £  H1+k(0)  and  w  £  H1~k{0)  =  4ffc_1(0)*.  We  can  chose  vn,wn  £  Hx(0), 
vn  — >  v  in  the  topology  of  H1+k(0)  and  wn  — >  w  in  the  topology  of  H1~k(0). 
Then 

B(wn,vn)  =  -  /  wnAvndx+  /  (wn\do)dvvndS(x). 

Jo  J  dO 

Next,  we  notice  that  Avn  Av  in  Hk~1(0),  and  hence  f0wnAvndx  — >  (w,  Av). 
Similarly,  d„vn  — >  dvv  in  Hk~1/2{dO).  Now,  if  the  sequence  B(wn,vn )  had  a  limit 
that  depended  only  on  w  and  v,  then  it  would  follow  that  the  sequence  of  traces 
Wnloo  would  have  a  limit  depending  only  on  w.  In  turn,  this  would  provide  by  a 
continuous  trace  map  H1~k(0)  — *  H1'2~k(dO),  k  >  1,  which  is  known  not  to  be 
possible. 

For  the  reasons  just  explained,  we  introduce,  following  an  idea  from  [22],  the 
space  H~s(0)  :=  H~s(0 )  ®  H~s~1^2(dO),  where  s  £  Z,  s  >  2.  Intuitively,  the 
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second  component  (  of  an  element  u  =  (it,  Q  in  H  s(0)  should  be  thought  of  as 
some  sort  of  trace  at  the  boundary  of  u.  We  then  define 

B  :  Hl~k(0)  x  H1+k(0)  ->  C  by 

(52)  B(u,V)  =  -(«,  AV)  +  <C,^F|  d0)do,  where 

u  =  («,  C)  e  H1~k(0)  ©  Hxl2~k(dO)  =  k£  N. 

With  this  definition,  we  can  now  introduce  weak  solutions  of  Equation  (1)  for 
r  =  1  —  k,  k  £  Z+. 

Definition  4.4.  Let  g  £  H~1/2~k(d 12),  k  £  Z+.  We  say  that  it  £  H1~k(fl)  satisfies 
the  Equation  (1)  (i.e.,  —A it  =  0,  dvu  =  g  £  H~1^2~k(dfl))  in  weak  sense  (or  that  u 
is  a  weak  solution  of  the  Equation  (1))  if,  and  only  if,  there  exists  (  £  H1^2~k(dfl) 
such  that 

B(u,  V)  =  ( g ,  V)ao,  where  u  =  (u,  Q 

for  all  V  £  H1+k(n). 

The  pair  u  =  (u,()  above  will  also  be  called  a  weak  solution  of  Equation  (1). 

Remark  4.5.  The  above  definition  of  weak  solutions  of  the  Neumann  problem 
generalizes  the  classical  definition.  Indeed,  if  g  is  regular  enough  so  that  u  £  H2(<D) 
is  a  classical  solution  of  the  Neumann  problem  with  data  g  (i.e.,  of  Equation  (1)) 
then 

B(u,V)  =  B(u,V)  =  (g,V)ao 

for  all  V  £  H1+k(tt),  where  u  :=  ( u ,  u|an).  Hence  u  is  a  weak  solution  of  Equation 
(1)  also  in  the  sense  of  Definition  4.4. 

We  shall  need  the  following  simple  observation. 

Lemma  4.6.  Let  u,u\  £  Hl~k(0 )  =  H1~k(0)  ©  H1/2~k(dO),  u  =  (u,  ()  and 
Ui  =  (iti,Ci)-  Then  B(u,V )  =  B(ii\,V)  for  all  V  £  H1+k(Q)  if,  and  only  if, 
u  —  Ui  =  c  and  (  —  (i  =  c,  where  c  is  a  constant. 

Proof.  The  subspace  V  :=  {(— AV,  dvV)}  C  Hk~1(0)  ©  Hk~1^2(dfl)  has  codimen¬ 
sion  one,  by  the  solvability  conditions  of  the  Neumann  problem.  Moreover,  the 
solvability  conditions  for  the  Neumann  problem  (recalled  in  the  proof  of  Proposi¬ 
tion  4.9)  show  that  the  annihilator  of  V  is  the  vector  (1,1).  Our  assumption  is 
equivalent  to  the  fact  that  w  :=  u\  —  u  satisfies  (w,u>i)  for  any  w\  £  V.  Since 
Hl~k(0)  is  the  dual  of  Lffe_1(0)  ©  Hk~1^2(dfl),  it  follows  that  u\  —  u  must  be  a 
multiple  of  the  non-zero  vector  (1, 1)  spanning  the  annihilator  of  V.  □ 

From  the  above  lemma,  Lemma  4.6,  we  obtain  the  following  corollary. 

Corollary  4.7.  Ifu  =  (it,  £)  is  a  solution  of  Equation  (1)  in  weak  sense  (Definition 
4-4)  then  f  is  uniquely  determined  by  u.  Moreover,  if  u\  =  (tti,£i)  is  another 
solution  of  Equation  (1),  then  u±  ~  u  is  a  constant. 

The  usual  properties  of  the  solutions  of  the  Neumann  problem  are  therefore 
satisfied.  The  above  corollary  also  allows  us  to  define  u\go  ■=  C  if  (UX)  is  a  weak 
solution  of  Equation  (1)  for  some  g  £  H~1!2~k(dO).  By  abuse  of  terminology,  we 
shall  say  that  u  £  H1~k(0 )  is  a  weak  solution  of  Equation  (1)  if  there  exists  (  £ 
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H  1/2  fc(i dO)  (uniquely  determined  by  the  above  Corollary)  such  that  u  =  (u,Q 
is  a  weak  solution  of  that  equation.  In  this  case,  we  shall  also  write 

(53)  B(u  +  Ul,V)  :=  B{u,V)  +  B{Ul,  V) 

for  any  u-\  G  if1(0).  This  is  needed  in  Equation  (54)  below.  See  also  (63). 

Let  us  G  S  be  the  GFEM-approximation  of  the  solution  of  Equation  (1)  (as 
defined  by  Lemma  4.1).  Also,  let  u  =  (m,  C)  be  a  weak  solution  of  Equation  (1). 
The  extension  (53)  of  the  definition  of  the  bilinear  form  B  is  useful  because,  for 
example,  it  allows  us  write  that 

(54)  B(u  —  us,  Vs)  =  0,  for  all  Vs  G  S. 

The  relation  (54)  is  then  equivalent  to  the  Definition  4.2. 

Let  X  and  Y  be  normed  spaces  with  norms  ||a:||x  and  ||y||r  and  let  B\  :  Xx  Y  — > 
C  be  a  bilinear  form.  Recall  that  B\  is  said  to  be  continuous  if,  and  only  if,  there 
exists  C  <  oo  such  that  \B\{x,y)\  <  C'||x||.x’||2/||y  f°r  all  x  G  X  and  y  G  Y.  We 
shall  need  the  following  result. 

Theorem  4.8.  Let  X  and  Y  be  reflexive  Banach  spaces  with  norms  ||x||x  and 
||j/||y.  Also,  let  B\  :  X  x  Y  — >  C  be  a  bilinear  form.  Assume  that 

(i)  B\  is  continuous; 

(ii)  There  exists  7  >  0  such  that 

inf  sup  \Bi(x,  y)\  >7; 

IMU=i  ||j/||y.<i 

(in)  sup  \Bi(x,  y)\  >  0  whenever  y  ^  0. 

IMIx<i 

Then  for  any  continuous  functional  F  :  Y  — ►  C  there  exists  a  unique  x  G  X  such 
that  F(y)  =  Bi(x,y),  for  all  y  G  Y.  Moreover,  we  have  ||x||  <  ||F||/7. 

This  theorem  is  a  generalization  of  the  well  known  Lax-Milgram  Lemma.  For  a 
proof,  see  Theorem  5.2.1,  page  112,  of  [1].  The  proof  in  that  book  is  an  adaptation 
of  the  proof  in  [17]  and  [19].  See  also  [18],  page  294. 

We  now  check  that  our  form  B  satisfies  the  conditions  of  Theorem  4.8. 

Proposition  4.9.  Let  X  C  Hl~k{0)  :=  141”fe(C>)  ©  Hl/2~k{dO)  consist  of  the 
pairs  u  =  (it,  £)  satisfying  (it,  1)  +  (£,  1)30  =  0.  Also,  let  Y  C  H1+k((D)  consist  of 
the  functions  V  such  that  ( V, ,  1)  :=  J0  V ( x)dx  =  0.  Then  the  restriction  of  the  form 
B  of  Equation  (52)  to  X  x  Y  satisfies  the  conditions  of  Theorem  f.8  for  k  >  1. 

Proof.  The  bilinear  form  B  is  immediately  seen  to  be  continuous  by  definition  of 
B  and  its  definition  (Equation  (52))  and  by  the  definition  of  our  negative  order 
Sobolev  spaces.  Therefore  Condition  (i)  in  Theorem  4.8  is  satisfied. 

To  check  that  Condition  (ii)  in  Theorem  4.8  is  satisfied,  we  shall  use  the  well 
posedness  of  the  Neumann  problem  on  O.  Namely,  we  shall  use  the  fact  that,  for 
any  v  G  Hk~l(0)  and  any  g  G  Hk~1^2(dO)  satisfying 

(l,v)  +  (l,g)g0  ■■=  /  v(x)dx+  /  g(x)dS(x)  =  0, 

J  n  Jon 

there  exists  a  unique  V  G  Hk+1(0)  such  that 

(55)  -AC  =  v,  dvV  =  g ,  and  (C,  1)  =  0. 
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Moreover,  there  exists  Cq  >  0  such  that 

(56)  II  ^11  Affc+i  (o)  <  C'si  ( 1 1  1 1  Hk-l(o)  +  ll'/llfffc-l/2(ao))- 

Let  now  u  =  (u,  £)  £  X  c  H1-k(0)®H1/2~k(d0)  =  (tf*-1^)©#*"1/2^©))*. 
By  definition,  there  exist  v  £  Hk~1(0)  and  g  £  Hk~1^2(0 ),  not  both  zero,  such 
that 

(57)  (u,{v,g)}  :=  (u,v)  +  ( C,g)ao  >  ^INI  (JMli^-qo)  +  llsll^-i /3(ao)  j  > 

where  ||u||2  :=  IM|^i_fc(0)  +  \\C\\%i/2-k^gc)y  and  the  last  paranthesis  stands  for  the 
norm  of  the  pair  (v,g)  in  the  Hilbert  space  Hk~1(0)  ©  Hk~1^2(dO). 

If  we  replace  v  with  v  +  X  and  g  with  g  +  A,  where  A  denotes  the  constant  function 
equal  to  A  £  C,  then  the  pairing  {u,v)  +  (( ',g)oo  does  not  change,  since  u  £  X. 
Moreover,  choosing  A  such  that  ||i>  +  A || ^fc_ x +  ||g  +  X\\2Hk_1/2^do^  is  minimal 
means  replacing  the  pair  ( v ,  g)  with  its  projection  onto  the  orthogonal  complement 
of  (1,1)  (we  have  used  here  Equation  (10)).  This  choice  will  not  affect  Equation 

(57) .  We  can  thus  assume  that  (v,g)  _L  (1,1).  Therefore  (l,i>)  +  (1  ,g)oo  =  0,  and 
hence  the  Neumann  problem  with  data  (v,  g)  is  solvable.  Let  us  chose  then  V  as  in 
Equation  (55). 

With  u  =  (u,  C)  €  X  C  H1~k(0)  ©  H1^2~k(dO),  as  above,  we  obtain 
B(u,  V )  =  -(u,  AH)  +  (C,  dvV\ do)oo  =  {u,  v)  +  (C,  g)ao 

>  2 INI  (IMI/ffc-1(e>)  +  \\9\\Hk-l/2(dO))  !  —  C\\V\\2Hk+l(oy 

This  verifies  Condition  (ii)  of  Theorem  4.8. 

Finally,  to  check  Condition  (iii)  of  Theorem  4.8,  let  V  £  Y  C  H1+k(0)  be  such 
that  B(u,  V)  =  0  for  all  u  £  X  C  fL1_fc(0).  Then  the  definition  of  the  space  X  as 
the  orthogonal  of  (1, 1)  shows  that  —AV  =  c  and  dvV  =  c,  for  some  c  £  C.  Green’s 
formula  gives 

0  =  -(1,  AH)  +  (1  ,dvV)ao  =  c(vol(H)  +  vol(cftl)), 

and  hence  c  =  0.  From  this  we  next  obtain  that  V  is  a  constant.  Since  V  £  Y,  this 
constant  must  also  be  zero.  □ 

We  therefore  obtain  from  Theorem  4.8  that  the  Neumann  problem,  Equation  (1) 
has  a  weak  solution  for  any  g  £  H~1//2~k(dO). 

Theorem  4.10.  Let  g  £  H~x/2~k{dO)  satisfy  (g,l)eo  =  0-  Then  there  exists 
u  =  (it,  C)  £  H1^k(0)  satisfying  the  Equation  (1)  in  weak  sense.  This  solution  is 
uniquely  determined  if  ( u ,  1)  =  0  and  then  it  satisfies 

llMlln1-fc(o)  +  IICI|j?i/2-»'(ao)  <  C'oll5lli4-i/2-f=(ao), 

for  a  constant  that  depends  only  on  O. 

Conversely,  let  u  =  (it,  £)  be  a  weak  solution  of  Equation  (1).  By  taking  Av  = 
0  and  v\go  arbitrary,  it  follows  from  the  definition  of  the  weak  solutions,  that 
llfl'll i/2— fe(0O)  —  C|ICIIffi/2-fc(ae>).  Then,  by  taking  v  =  0  on  dO  but  Av  arbitrary, 
we  again  obtain  from  the  definition  that 

(58)  IICIIffi/2-fc(0e>)  <  C\\u\\Hi-k(0y 
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Finally,  this  gives  that 

(59)  llsfllff-i/2-fc(3e>)  -  C|lull m-k(0), 

as  in  the  case  k  <  0  (note  however  that  in  the  case  k  >  0,  we  also  have  A u  =  0  in 
distribution  sense) . 

5.  Approximate  solution  of  the  Laplace  equation  with  distribution 

BOUNDARY  CONDITIONS  USING  THE  GFEM 

We  shall  consider  the  same  setting  as  in  the  previous  sections.  For  instance, 
S  =  Sv,  v  £  Z+,  will  be  the  GFEM-space  associated  to  any  of  the  data  T,v 
satisfying  the  assumptions  of  Subsection  2.3.  Also,  ft  =  O,  a  smooth,  bounded 
domain,  as  in  the  previous  section. 

From  now  on,  we  shall  fix  the  weak  solution  u  £  H1~k(0),  k  £  N,  k  <  m  —  1,  of 
the  Equation  (1)  satisfying  (u,  1)  =  0.  Recall  from  Definition  4.4  that  this  means 
that  there  exists  £  £  H1/2~k{dO)  such  that  u  =  (u,  £)  is  a  weak  solution  of  the 
Equation  (1)  in  the  sense  that 

(60)  B(u,  V)  =  (; g ,  V)d0,  for  all  V  £  H1+k(0), 

where  B{u,V)  =  —(u,AV)  +  (£,  dvV\go)dOi  see  Equation  (52).  Here  the  data 
g  £  H~l/2~k{0)  is  also  fixed. 

With  u,  £,  and  g  as  in  the  paragraph  above,  we  shall  define 

(61)  u\dn  :=  £  and  dvu\au  ~  9- 

We  shall  think  of  dvu  as  the  normal  derivative  of  u  in  the  direction  of  the  outer 
normal  at  the  boundary.  Recall  from  Equations  (58)  and  (59),  that  £  and  g  above 
depend  continuously  on  u.  In  the  forthcoming  paper  [6],  we  shall  compare  our 
definition  of  a  solution  of  Equation  (1)  with  the  definitions  in  [15,  26]  or  [28].) 

We  shall  also  assume  that 

(62)  (u,  1)  =  (g,  l)do  =  0, 
which  guarantees  the  existence  of  u  and  that 

llull Hi~k(o)  +  ||Cllui/2-«'(ao)  <  C'o\\g\\H~i/2-k(dO)^ 

by  Theorem  4.10. 

These  condition  of  Equation  (60)  imply  that  (u,  A(j>)  =  0  for  all  £  C%°(0),  that 
is,  that  A u  =  0  in  the  sense  of  distributions  on  O.  Since  £  is  determined  by  it,  we 
can  write  B{u,v)  :=  B{u,v).  More  generally,  we  shall  write 

(63)  B(u  +  Hi,  v)  =  B(u,  v)  +  B(ui,  v), 
whenever  u\  £  I/1(G).  See  also  Equation  (53).  In  particular, 

B(u  —  uv,  v)  =  0, 

where  u„  =  us„  £  S „  is  the  GFEM-approximation  of  u,  see  Definition  4.2  and 
Equation  54.  This  is  in  agreement  with  the  results  in  [15]  (especially  Theorem  6.5) 
on  traces  of  functions  w  such  that  Aw  is  regular  enough. 

We  have  the  following  estimate. 

Lemma  5.1.  With  u  as  in  Equation  (60)  above,  we  have 

\B(u,v)\  :=  \B(u,v)\  <  C\\u\\Hi-k{0)\\v\\Hi+k(0)l 

for  any  v  £  H1+k(0)  and  a  constant  C  depending  only  on  O. 
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Proof.  By  definition,  using  also  Equation  (58),  we  have 

(64)  \B(u,v)\  =  |  -  (u,  Aw)  +  {u\ao,dvv)do\  <  |(«,  Au)|  +  \{u\ao,  dvv)ao\ 

<  |M|ffi-fc(o)l|Au||ff-i+fe(o)  +  \\u\\Hi/2-k(d0')\\duv\\H-i/2+k(d0) 

<  C\\u\\Hl-k(0)\\v\\Hl+k(0y 

This  completes  the  proof.  □ 

We  continue  with  more  lemmas.  We  have  the  following  “inverse  property.” 
Lemma  5.2.  We  have  that 

(65)  \\w\\ h*(0)  <  -#«(£>), 

for  0  <  s  <  t  <  m  and  w  £  Sv,  for  a  constant  B  independent  of  v. 

Proof.  For  s  =  t,  the  result  is  tautologically  true  with  B  =  1.  For  s  =  0,  the  result 
is  given  by  Proposition  2.5.  Since  O  is  smooth,  the  general  case  by  interpolation 
using  the  results  of  [15]  on  the  interpolation  properties  of  Sobolev  spaces  on  smooth, 
bounded  domains.  □ 

In  the  proofs  below,  we  shall  occasionally  denote  h  =  hu. 

Lemma  5.3.  The  GFEM-approximations  uu  satisfy 

IKHffbO)  <  Chv  k\\u\\Hx-k(0) 

for  a  constant  C  depending  only  on  O  (so  C  is  independent  of  v). 

Proof.  The  Poincare-Friedrichs  inequality,  Lemma  5.1,  and  Lemma  5.2  give 

ll^lljrqo)  —  C  B(uu,  u„)  =  C B(u,  11^)  <  1 1  ^  1 1  x — fc  (c>)  1 1 1 1  i+fc  (O) 

<  Ch  k\\u\\Hi-k(C,')\\uv\\Hi(o), 

where  in  the  last  inequality  we  have  used  Lemma  5.2  for  s  =  1  and  h  =  □ 

This  gives  the  following  corollaries. 

Corollary  5.4.  Let  k  <  A.  Then  the  GFEM-approximations  u„  satisfy 

IWI H1~k(0)  <  C\\u\\Hi-k(Q) 

for  a  constant  C  depending  only  on  O.  In  particular,  uu  depends  continuously  on 
u  and  ||u  -  u„||iri-*(0)  <  C\\u\\Hi-k^0y 

Proof.  The  result  is  well  known  for  k  =  0  since  u v  is  the  B-orthogonal  projection 
of  u  onto  S  (this  is  Cea’s  Lemma,  see  [9,  10]).  We  shall  therefore  assume  that  k  >  1. 
Let  h  =  hu. 

Let  v  £  Fffc_1((D)  be  arbitrary.  Let  c  £  C  be  such  that  f0(v  —  c)dx  =  0.  Then 
we  can  find  V  £  Hk+1(0)  such  that 

-AV  =  v-c,  f  Vdx  =  0,  dvV  =  0,  and  ||V||#fc+i(e,)  <  C\\v\\Hk-it0\, 

Jo 

where  C  is  a  constant  depending  only  on  O.  Also,  chose  w  £  S  such  that 
<  Clip'll .F/fc+qo)  and  11^  —  w\\Hi(o)  <  Chk\\V\\Hk+i(0y 


DISTRIBUTIONS 


27 


This  is  possible  by  Theorem  3.2.  Then 

(uv,  v)  =  (uv,  v  -  c)  =  -(u„,  AV )  =  -(u„,  AV)  +  (uv\do,  dvV)ao  =  B{uv,  V) 

=  B(uv,  w)  +  B(uv,  V  —  w)  =  B(u,  w)  +  B(u„,  V  —  w). 
Using  also  Lemmas  5.1  and  5.3,  this  gives 

\(uu,v}\  <  C'||w||fl-i-fc(n)||ui||ffi+fc(n)  +  ||w1/||ffi(n)||V'  -  w||hi(q) 

<  C,||w||_H-l-fc(n)||1l/r||i7fc+l(n)  +  Ch  k\\u\\Hi-k^Q-)hk\\V\\Hk4ri^') 

<  C,|lullff1-fc(n)l|1l/lli7fc+1(n)  —  C|Mli?1-fc(n)IMIirfc-1(n)- 

This  gives  the  result  since  ||uj/||ifi-fc(r2)  :=  sup  \{uv,  u)|/|M|i?fc-i(n),  v  ^  0.  □ 

Similarly, 

Corollary  5.5.  We  have  \\uu\do\\H1/2-k(dO)  —  ^\\u\\H1~k(0)  for  a  constant  C 
depending  only  on  O .  In  particular,  ||(w  —  uv)\do^H1/*~k{dO)  —  C|lullff1-'E(0) • 

Proof.  The  proof  is  similar  to  that  of  the  previous  corollary.  Let  v  £  H~1^2+k(dO) 
be  arbitrary.  Let  c  £  C  be  a  constant  such  that  Jgo  vdS  =  J0  cdx.  Then  we  can 
find  a  unique  W  £  H1+k(0)  satisfying 

AW  =  c,  [  Wdx  =  0,  dvW  =  v,  and  ||W||hW(q)  <  C\\v\\H-i/2+k,d0), 

Jo 

for  a  constant  C  >  0  depending  only  on  O.  Using  also  Theorem  3.2,  we  choose  w  £ 
S  such  that  <  C\\W\\Hk+i^  and  \\W  -  w\\Hi^  <  Chk\\W\\Hk+i{Qy 

Then,  using  also  {u„,  1)  =  0,  we  obtain 

{uv\ dO,v)dO  =  (uv\oo,dvW)do  =  {uv,  AW)  +  B(uv,  W)  =  B(uv,  W) 

=  B{uv ,  w)  +  B{uv ,  W  —  w)  =  B(u,  w)  +  B(u„,  W  —  w). 
Using  Lemmas  5.1  and  5.3,  we  then  obtain 

\{uu\oo,v)do\  <  U||u||Hi-fc(n)||ui||ffi+fc(n)  +  ||u,/||fl-i(n)||W  -  HI-tfRn) 

<  C|lull77l-fc(n) l|W/II.F/l+fc(n)  +  Ch  k\\u\\Hi-k^')hk\\W\\Hi+k^cl) 

<  C|lull/f1-fe(fi)l|Ur||ff1+fe(e>)  <  C\\u\\Hi-k(Qj-)\\v\\H-i/2+k(d0), 
which  completes  the  proof  in  view  of  the  definition  of  ||w1/||  H1/2~k(Qy  D 

We  now  establish  our  main  approximation  results.  Recall  that  u  £  H1~k(0)  is 
the  weak  solution  of  A u  =  0,  dvu  =  g  £  H~1'2~k{dO)  (Equation  (60))  satisfying 
(u,  1)  =  0  and  that  uv ,  v  £  Z+,  are  the  GFEM-approximations  of  u,  where  S  = 
Sv  are  the  GFEM-finite  element  spaces  associated  to  a  sequence  of  data  E„  = 
{ivf,  ff,  tyj,  x>*"}  satisfying  the  assumptions  of  Subsection  2.3.  The  constants  C  >  0 
below  are  assumed  to  be  independent  of  Sv. 

Proposition  5.6.  Assume  that  fc  +  7  <  A  and  k  <m  —  1,  7  £  Z+.  Then  the  error 
u  —  uu  £  H1^k(0)  satisfies 

llu  —  wll/fi-fc— qn)  <  C,^llulli71-fc(n)) 

with  a  constant  C  independent  of  v. 
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Proof.  Let  v  €  H  1+k+1(0)  be  arbitrary.  Let  c  be  a  constant  such  that  (v  —  c,  1)  = 
0.  Then  there  exists  a  unique  V  £  H1+k+1  (O)  such  that 


for  a  constant  C  >  0  depending  only  on  O. 

Then,  for  any  w  £  S, 

(66)  (u  —  uv,  v)  =  (u  —  u„,  v  —  c)  =  ~(u  —  uv,  Ay) 

=  ~(u  -  uv,  Ay)  +  ((u  -  uv)\dO,dvV)do  =  B(u  -  uv,  V )  =  B(u  —  uv,V  —  w) 

=  ~(u  -  Uv,  A(y  -  w))  +  ((U  -  Uv)\dO,9vW)dO- 

Since  fc  +  7  <  A  and  k  <  m—  1,  the  assumptions  of  Theorem  3.2  are  satisfied,  so  we 
can  chose  w  £  S  such  that  ||io||#i+k+7(f2)  <  C||y||jji+fc+7(f2)  and  \\V  —  ^ 

Ch't\\V\\Hi+k+-,{n).  In  particular, 

\\dvw\\H~i/2 +k(ao)  =  \\dv(V  -  w)\\H-i/2+Hd0)  <  /?7||y||/i-i+fc+7(o). 

From  || y || i/i+fc+7 (n)  <  C'||i'||H-1+fc+^(0))  Corollaries  5.4  and  5.5,  and  Equation  (66), 
we  then  obtain, 

\(u-Uv,v)\  <  \\u-Uv\\Hi-k(Q.)\\V —  w\\Hi+k(a.)  +  \\u  —  uA\H1-k(Q.)\\dvW\\H-i/2+k(dO) 

<  C'ft7||u||/i-i-t(o)||u||fl--i+fc+7(e)). 

The  proof  is  complete.  □ 

We  then  obtain. 

Proposition  5.7.  Assume  that  the  local  approximation  spaces  'ly  contain  the 
polynomials  of  degree  A  >  k  +  'y,  7  £  Z+,  and  let  A  <g  B  <g  O  be  open  subsets.  Then 
for  hv  small  enough  and  k  +  7  <  m  +  1,  k  <  m  —  1,  and  A  >  l  >  1,  we  have 

!|w  -  uv\\m(A)  <  Chlv\\u\\Hi+x(B)  +  CKl\\u  -  Uv\\h%^{b), 

where  the  constant  C  depends  on  the  structural  constants  and  the  distance  from  dA 
to  dB,  but  not  on  v. 

Proof.  Assume  first  that  u  £  B1(0)  and  let  us  replace  m  with  fc+7  — 1  in  Theorem 

3.12,  which  is  possible  since  fc+7  <  m+1  and  we  can  always  decrease  in  in  Theorem 

3.12.  The  use  Theorem  3.12  for  A  <s  Bi  C  B,  for  some  admissible  open  subset  Bi 
such  that  B[  C  B.  The  result  then  follows  from  Proposition  5.6  and  from 

inf  ||U  —  x\\h1(b1)  <  Chl\\u\\Ht+i(B[)  <  C'/i,||u||ffi+i(B). 

In  general,  for  u  £  Pl1~k{0),  the  result  follows  by  continuity,  using  Corollary  5.4 
and  Equation  (3).  □ 

We  keep  all  our  previous  assumptions,  the  most  important  of  which  were  recalled 
before  Proposition  5.6.  In  particular,  the  assumptions  of  Subsection  2.3  remain 
valid.  Also,  recall  that  the  local  approximation  spaces  \Ey  are  assumed  to  contain 
the  polynomials  of  degree  A.  By  taking  l  =  7  and  using  again  Equation  (3),  we 
obtain  the  following  error  estimate  for  the  GFEM-approximation  u„. 


DISTRIBUTIONS 


29 


Theorem  5.8.  Let  u  £  H1^k(0)  be  a  weak  solution  of  the  Neumann  problem 
A u  =  0,  duu  =  g  £  H~1'2~k(dO)  and  £  Sv  be  its  GFEM-approximation,  where 
the  sequence  Sv  satisfies  the  assumptions  of  Subsection  2.3  with  k  <  min{A,m—  1}. 
Also,  let  Aq  <e  O  be  an  open  subset  of  the  smooth,  bounded  domain  O.  Then  for  hv 
small  enough. 

||w  -  uv\\Hi(Ao)  <  Chl\\u\\Hi-k(0),  7  =  min{A, to  +  1}  -  k. 

The  constant  C  above  depends  only  the  structural  constants  and  the  distance  from 
OAo  to  dO,  but  not  on  v. 

Remark  5.9.  The  constant  C  >  0  above  may  depend  on  7  and  on  the  structural 
constants  A,  Cj,  k,  a,  A,  and  m,  as  well  as  on  Aq  and  O. 

In  particular,  since  ||w||ffi-<=(0)  <  C\\g\\H-i/2-k(do)i  by  Theorem  4.10,  we  obtain 
the  following  corollary. 

Corollary  5.10.  Under  the  assumptions  of  Theorem  5.8,  we  have 
\\u  -  uu\\Hi(Ao)  <  Ch1\\g\\H-i/2-k(d0), 

with  C  >  0  independent  of  v  and  /i„  small  enough. 

Remark  5.11.  The  condition  that  hu  be  small  enough  is  not  an  essential  restric¬ 
tion.  Indeed,  by  increasing  the  constants  C  in  the  above  results,  we  can  drop  the 
requirement  that  hv  be  small  enough. 

6.  Polynomial  local  approximation  spaces 

In  this  section  we  shall  verify  that  the  Conditions  A-D  are  verified  if  we  choose 
'Pj  =  Qm,  to  be  the  space  of  polynomials  of  degree  <  m,  1  <  to,  and  the  boundary 
of  fl  is  piecewise  smooth,  for  a  suitable  covering  and  subordinated  partition  of  unity. 
Some  results  in  this  section  are  either  elementary  or  well  known.  We  include  them 
nevertheless  for  the  benefit  of  the  reader  and  for  completeness. 

For  any  ball  B  of  radius  r,  we  shall  denote  by  tB  the  ball  with  the  same  center 
as  B  and  radius  tr. 

Lemma  6.1.  There  exists  a  constant  C  >  0,  depending  only  on  n,  to,  and  M, 
such  that  for  any  ball  B  C  Mn,  for  any  Q  £  Qm,  and  for  any  t  £  (0 ,M],  we  have 
WQWLHtB)  <  C\\Q\\L2(B). 

Proof.  For  any  fixed  B,  Q  1— >  || <511  z,2 (*s)  and  Q  |— >  ||Q||l2(b)  are  two  norms  on 
the  finite  dimensional  space  Qm  of  polynomials  of  degree  <  to,  and  hence  they  are 
equivalent.  This  gives  the  result,  except  the  independence  of  C  on  B  and  t.  But 
all  balls  are  affine  equivalent  and  the  L2-norm  is  scaled  by  the  (square  root  of  the) 
determinant  of  the  matrix  of  the  affine  transformation.  Thus  the  constant  C  can 
be  chosen  to  be  the  same  for  all  balls  B.  □ 

This  gives  immediately  the  following  corollary. 

Corollary  6.2.  There  exists  a  constant  C  >  0,  depending  only  on  n,  in,  and  M 
such  that  for  any  ball  B  C  R”,  for  any  polynomial  Q  €  Qm >  and  for  any  t  £  (0,  M] , 
we  have  \Q\Hi(tB)  <  C\Q\hi(B)  and  \\Q\\Hi(tB)  <  C\\Q\\hi(B),  0  <1  <m. 

Proof.  Use  Lemma  6.1  for  all  derivatives  Q^a\  where  |a|  <1.  □ 

We  now  establish  to  the  following  “inverse  property.” 
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Lemma  6.3.  There  exists  a  constant  C  >  0,  depending  only  on  n,  m,  and  a,  such 
that  WQ^Wl2^)  <  Crl-\a\\\Q\\Hi(B)  for  any  l  <  |a|  <  m,  any  Q  G  Qm,  and  any 
ball  B  of  radius  r. 

Proof.  Let  us  prove  first  the  result  for  1  =  0.  That  is,  we  need  to  prove  that 
\\Qia)\\L*(B)  <Cr-W\\Q\\L2{B). 

Let  Bi  =  Bi(0)  be  the  unit  ball  centered  at  0.  Then  Q  i— >  ||Q^||B2(Bl)  is  a 
semi-norm  on  Qrn ,  the  space  of  polynomials  of  degree  at  most  m,  and  hence  it  is 
bounded  by  the  norm  Q  i-»  ||<3||B2(Bl).  Thus  ||Q(a) Wl^b^  <  Ci\\Q\\L2(Bly  Let  L 
be  an  affine  transformation  mapping  B\  onto  the  ball  B  of  radius  r  consisting  of 
the  composition  of  a  translation  and  a  dilation  of  ratio  r  (so  det(L)  =  rn).  Then 

l|Q(Q)l|zG(B)  =  det(L)1/2||Q^  o  L||L2(Bl)  =  det(L)1/2r"H ||(Q  o  L)<“> ||B2(Bl) 

<  det(L)1/2r-l“l \\Q  o  L\\l2(Bi)  =  Cr~^  \\Q\\l^b ), 

for  any  Q  G  Qm. 

Assume  now  that  |a|  >  l  >  0.  Choose  (3  <  a,  \(3\  =  l.  Then 

||£>“-/3£>/JQ||i2(B)  <  (7r-l“-dl ||D/3Q||L2(b)  <  Crl-W\\Q\\Hl(B). 

This  completes  the  proof.  □ 

The  relevant  “inverse  property,”  implying  also  Condition  C,  now  follows. 

Proposition  6.4.  There  exists  a  constant  C  >  0,  depending  only  on  n,  m,  a,  and 

а,  such  that  ||Q(“)||B2(0)  <  Crl^a^\\Q\\Hi  ( B)  for  any  l  <  |a|  <  m,  any  Q  G  Qm, 
any  ball  B  of  radius  r,  and  any  set  O  contained  in  a~1B  and  star-shaped  with 
respect  to  B. 

Proof.  This  follows  from  Corollary  6.2  and  Lemma  6.3.  □ 

We  now  prove  the  following  elementary  lemma. 

Lemma  6.5.  Assume  the  data  defining  S  is  fixed.  Then  for  any  k  G  Z+  and  any 
open  sets  Aq  <g  A C  f l,  we  can  construct  admissible  open  sets  B\, . . . ,  Bk-i  such 
that  A0  =:  B0  <e  Bi  <e  B2  <s  ■  ■  ■  <e  Bk-i  <e  Bk  :=  Aj~  and  Cdist(dBj,dBj_i )  > 
9/k,  where  9  :=  dist(dA0,dAfc),  provided  that  Ch  <  9,  where  C  depends  on  k  and 
n  only  (in  particular,  C  is  independent  of  h). 

Proof.  Take  k  =  2,  for  simplicity.  The  general  result  is  proved  similarly  or  by 
iterating  this  case.  Let  (7  =  4.  Let  U  be  the  union  of  all  open  sets  u>j  at  distance 
at  most  9/ 4  from  Ao.  Let  J  be  the  set  of  indices  j  such  that  4>j  /  0  on  U  and  let 
G  be  the  set  where  ‘Aj  =  1-  We  then  define  B i  to  be  the  interior  of  G.  □ 

б. 1.  Partition  of  Unity.  We  show  in  this  subsection  that,  for  a  piecewise  C 1  do¬ 
main  ft,  we  can  choose  a  sequence  Y,v  =  {ujf,  (j)1',  T'' ,  v  G  Z+,  of  GFEM- 

data,  U jixf  =  ft,  satisfying  the  assumptions  of  Subsection  2.3.  Other  examples  of 
sequences  Ytv,  possibly  better  suited  for  numerical  implementation,  will  be  included 
in  the  forthcoming  paper  [7],  where  we  will  also  discuss  the  numerical  implementa¬ 
tion  of  the  GFEM  for  boundary  value  problems  with  distributional  data. 

Theorem  6.6.  Let  O  C  R"  be  a  bounded  open  subset  with  piecewise  Cl -boundary 
and  all  angles  >  0.  Then  there  exist  structural  constants  A,  Cj,  a,  k,  A,  and  m 
such  that,  for  any  small  enough  h  >  0,  we  can  construct  a  partition  of  unity  {4>j} 
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subordinated  to  the  covering  {loj}  and  satisfying  all  the  assumptions  of  subsection 
(2.3),  in  particular,  satisfying  conditions  A(h),  B,  C,  and  D. 

Proof.  Let  us  first  triangulate  12  with  (possibly  curvilinear)  triangles  Tj  with  di¬ 
ameters  <  h/ 2  and  >  C\h  and  satisfying  the  angle  condition.  (Recall  that  this 
means  that  there  exists  an  angle  9  >  0  such  that  all  angles  of  the  triangles  are  >  9, 
with  9  independent  of  h.)  Only  the  sides  on  the  boundary  of  12  are  allowed  to  be 
curvilinear.  The  interior  sides  are  required  to  be  straight  segments.  This  is  possible 
since  12  has  a  piecewise  C 1  boundary  and  all  angles  are  >  0. 

We  can  then  chose  a  >  0  small  enough  and  independent  of  h  such  that  each 
triangle  will  contain  a  ball  of  radius  >  2ah.  Let  Xj  be  the  centers  of  these  balls 
and  let  to*  be  the  ball  with  radius  ah  and  center  Xj.  Then  let  u>j  =  B(xj ,  h)  fi  12, 
where  B(x,r )  denotes  the  open  ball  with  center  x  and  radius  r.  Our  construction 
shows  that  the  balls  2 w*  =  B(xj,2ah)  are  disjoint. 

Then  12  =  U f—i^Jj-  Let  <j>j(x)  =  r](\x—Xj\/h),  where  77  is  a  fixed,  smooth  function 
77  :  [0,oo)  — »  [0,1]  such  that  77(f)  =  1  for  t  <  1/2  and  77(f)  =  0  for  f  >  1.  Then 
W^Wl  00 (fi)  <  Cfc/(diamwj)fc,  for  k  =  |a|  <  m,  and  any  j  =  1, ...  ,7V.  Moreover, 
<f>j  =  1  on  Uj  and  hence 

E^1- 

Let  now  ipj(x)  =  77(1  a;  —  Xj\/(ah)),  so  that  ipj  has  support  in  2 lo*  and  is  equal 
to  1  on  u*.  Define  'ipj  =  1  —  ipj  and  let  ip  =  Ylf=1ipj-  Then  ip  =  0  on  Uw* 
and  || Oa(pj || £00 (fi(  <  Ck/(dia,m.coj)k,  for  k  =  |a|  <  m,  and  any  j  —  Let 

<pj  :=  ipcpj  +  ipj.  An  easy  verification  shows  that  we  still  have  JT  <pj  >  1,  because 

(i)  the  balls  2u>*  are  disjoint; 

(ii)  <pj  =  <pj  outside  U2w*, 

(iii)  <pj  =  1  on  2 lo*,  and 

(iv)  (pj  >  0  everywhere. 

Finally,  we  consider  the  Shephard  functions 


k 


Our  construction  shows  that  (pj  is  a  (k,  Cq,  C\, . . . ,  Cm)  partition  of  unity.  More¬ 
over,  (pj  =  1  on  u>j,  because  (pk(x)  =  6k j  for  x  €  uj*  (recall  that  6k j  =  1  if  k  =  j 
and  6k j  =  0  otherwise).  □ 

Some  assumptions  on  the  domain  12  in  the  above  theorem  are  necessary,  as  shown 
by  the  following  remark. 

Remark  6.7.  The  non-Lipschitz  domain 

12c  :=  {(x,y),  -x2  <  y  <  x2,  x2  +  y2  <  1,  x  >  0} 
will  have  no  covering  {u>j}  satisfying  the  Condition  A. 
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